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ABSTRACT 

SN 2004et is one of the nearest and best-observed Type IIP supernovae, with a progenitor detection as well as good photometric 
and spectroscopic observational coverage well into the nebular phase. Based on nucleosynthesis from stellar evolution/explosion 
models we apply spectral modeling to analyze its 140 - 700 day evolution from ultraviolet to mid-infrared. We find a Mzams = 15 
Mq progenitor star (with an oxygen mass of 0.8 Mq) to satisfactorily reproduce [O i] AA6300, 6364 and other emission lines of carbon, 
sodium, magnesium, and silicon, while 12 Mq and 19 M© models under- and overproduce most of these lines, respectively. This result 
is in fair agreement with the mass derived from the progenitor detection, but in disagreement with hydrodynamical modeling of the 
early-time light curve. From modeling of the mid-infrared iron-group emission lines, we determine the density of the "Ni-bubble" to 
p(t) 7 X lO"''' X (t/100 d)"^ g cm"', corresponding to a filling factor of / = 0.15 in the metal core region (V = 1800 kms"'). We 
also confirm that silicate dust, CO, and SiO emission are all present in the spectra. 
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. 1. Introduction 

. Stars with Zero Age Main Sequence (ZAMS) mass greater than 
• about 8 Mq end their lives as core-collapse supernovae (SNe). 
[ Over half of these even ts (per unit volume) are classified as 
. Type IIP dLi et al.ll201 ih . showing hydrogen lines, as well as a 
'3-4 month plateau in the light curve, implying the presence of 
[ a massive hydrogen envelope. After the plateau phase, the core 

■ of the SN becomes visible, glowing from radioactive input by 

■ ^^Co. As long as the ejecta remain opaque to the gamma-rays 
, emitted in the decay, the light curve follows the exponential de- 

■ cay of ^^Co, with an e-folding time of 1 1 1 .4 days. The spectrum 
[ evolves from being dominated by absorption lines and lines ex- 
. hibiting P-Cygni profiles superimposed upon a blackbody-like 

■ continuum to strong emission lines and a weaker continuum. As 
I the temperature falls, thermal emission eventually shifts into the 

■ infrared, but ultraviolet/optical features remain, caused by non- 

■ thermal ionizations and excitations. 

From hydrodynamical modeling, the proge nitors of Type IIP 
SNe are believed to be red supergiants (RSGs) (IChevalieilll976t 
lFalk& Arnet3ll977h . The rapidly growing image archives have 
improved the prospects of confirming this by identifying SN 
progenitors in imaging surveys. Such identifications have now 
been made for the Type IIP SNe SN 2003gd, SN 2004A, SN 
2004et SN 2005cs, SN 2008bk, SN 2008cn, and SN 2009md 
(ISmartlll2009l: lEhas-Rosa etaDl2009l: iFraseret alJl201lL and ref- 
erences therein). For SN 2003gd, the progenitor identification 
was verified by follow-up observations showing that the can- 
didate had disappeared (Mauiid & Smartt 2009)- In the cases 
where multiple-band detections were made, the progenitors have 
shown properties consistent with the RSG hypothesis, except for 



SN 20 08cn, where the progenitor was yellowQ di ias-Rosa et al.l 
I2OO9I) . 

The estimated luminosity of the progenitor allows its ZAMS 
mass to be determined from stellar evolution models. The 
resulting values are consistently found to be below 20 Mq, 
with a statistical analysis yielding a range of 8-5^] ^ - 16.5;^['^ 
Mq for the progenitor population dSmartt et al.ll2009l) . although 
dust effects may produce a s omewhat higher upper boundary 
(IWalm swell & Eldridgel20T2l) . The fate of stars in the 20-25 /30 
Mq range then remains unclear, as standard stellar evolution 
models predict also these to evolve to RSGs and explode as Type 
IIP SNe (e.g. iHeger & Lange^l2000l: iMevnet & Maede^l2003l: 
lEldridge et al.ll2008h . Models including rotation show that stars 
with initial mass greater than 20 Mq and time-averaged equa- 
torial velocities of >200km s"' could lose much of their hy- 
drogen envelope and mo ve bluewards in the HR diagram (e.g. 
iMevnet & Maedeil2003h . While this is immediately attractive as 
an explanation for the lack of luminous RSGs as Type IIP pro- 
genitors, one must accou nt for the broad obse rved distribution of 
equatorial velocities (e.g. lHunter et al.l[2008h as well as the lack 
of detection o f blue, luminous progenitors for Type Iln/Ib/Ic SNe 
(ISmarttll2009t) . 

Another method to derive information about the exploded 
star is through radia tion-hydrodynamical rn o deling of the bolo- 
metric light curve. iLitvinova & NadezhinI { 1981 [1981 pre- 
sented scaling relations from fits to a grid of Type IIP SN models 
allowing for the determination of the explosion energy, ejecta 
mass, and progenitor radius from the observed plateau length, 

' There is, however, the possibility that the detected source is a blend 
of two or more stars, and SN 2008cn differs from normal Type IIP SNe 
in its plateau length. 
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V-band magnitude, and photospheric velocity. However, the ne- 
glect of ^^Ni in these models renders them viable only for SNe 
with very small ^^Ni masses. [Eastman et al. ( 1994) found that 
a (t ypical) ^^Ni mass of 0.06 Mq extends the plateau by 40%, 
and iKasen & WooslevI (l2009l) found a prolongation by 20-30% 
for ^^Ni-masses of 0.05 - 0.1 Mq. Since the ejecta mass has 
a s trong scalin g with the plateau length (to the power 2.9) in 
the lLitvinova & Nadezhin ( 1985) relations, a significant overes- 
timate of the ejecta mass results. 

More reliable results can be obtained by fitting hydrody- 
namical models including ^*'Ni to the whole light curve evo- 
lution. Such modelin g has been undertaken for SN 19 99em 
dBaklanov et al. 2005; [Utrobinl 120071: i Bersten et al.l 1201 Ih. SN 
2003Z (lUtrobin et all l2007h. SN 2004et (lUtrobin & Chugail 
I2009h. SN 2005cs (IUtrobin& Chugail 120081) . and SN 2009kf 
(lUtrobin et alj20l"o|). obtaining ejecta masse s Mej - 14-28Mo. 
As dis cussed by Utrobin & Chugail (l2009l) and iMaguire et alj 
(1201 Ol MIO hereafter), these ejecta masses are generally too 
high to be consistent with the initial masses determined from 
direct observations of SN progenitors, as well as f or what is ex- 
pected fro m stellar evolu t ion in general. Recently, llns erra et alJ 
(1201 lb and llnserra et al. ( 20121). us i ng the radiation hydrody- 
namics code of Pumo & Zampieril (1201 lb . determined lower 



ing produced by the gamma-rays and positrons, which is the 
first step in the nebular phase modeling. The evoluti on of in- 
dividual lines were analyzed in a s e ries o f papers by IXu et al 



ejecta masses of Mej = 5 - 7.5 and Mej = 8-12 Mq for the 
Type IIP SNe 2007od and SN 2009bw. 

A third method for diagnosing the progenitor mass is through 
late-time spectral modeling. In the nebular phase (f > 150 days), 
the inner ejecta become visible, and the various nuclear burning 
zones can be analysed. Stellar evolution models predict the metal 
core mass to s trongly increase with progenitor ZAMS mass (e.g. 
IWooslev & Weaver. 1995l) , making it possible to distinguish be- 
tween different progenitors by determining the nucleosynthesis 
yields. Emission lines of carbon, oxygen, neon, sodium, magne- 
sium, silicon, and sulphur are the main signatures that can po- 
tentially constrain the progenitor mass. 

Despite being the most common core-collapse SN class, 
there are so far only a handful of Type IIP SNe that have exten- 
sive temporal and spectral coverage in the nebular phase, and not 
much information has as yet been derived from these observa- 
tions. The best existing datase t s are for SN 1990E (ISchmidt et"aD 
[l99l iGomez & LopezllT99l iBene tti et al][T994l) . SN 1999em 
(iLeonard et all l2002at lElmhamdi et all l2003l). SN 1999gi 



Leonard et al.' '2002b), SN 2002hh ("Pozzo et al.' '2006), SN 
2004dj (Chugai et al. 2005), SN 2 004et (Sahu et al. 20 06, MIO) 
SN 2 005CS (e.g. iPastorello et alJl2009l). 2007it (lAndrews et alJ 
120111) . and 2007od (Jnserra et alJ 1201 ll) . Of these, only SN 
1999em was well monitored also in the near-infrared (NIR). 

Apart from the sparsity of data, spectral analysis needs 
complex tools to model the transformation of energy from ra- 
dioactive decay to emergent UVOIR radiation. Temperature and 
NLTE solutions for a large number of zones and elements, 
including non-thermal heating/ionization/excitation rates, and 
multi-line radiative transfer are the fundamental ingredients that 
must be calculated, together making up millions of constraints 
whose self-consistent solution must be iteratively sought. In ad- 
dition, molecules, dust, and time-dependent effects may cause 
further complications. 

So far, SN 1987 A is the only Type II SN for which de- 
tailed spectral modeling has been undertaken. Despite being 
a Type Ilpec rather than a Type IIP, it was probably similar 
to a Type IIP in the nebular phase, since the nucleosynthe- 
sis is l argely unaffected by th e lat e-time evolution of t he en- 
velope. IXu & McCravl ( Il991h and iKozma & FranssonI d 19921) 
obtained the solutions for the ionization, excitation, and heat- 



('I99j (hydrogen), iLi & McCravl d 19921) ( oxyge n), iLi et al 
il99% (iron, coba l t, nick el), iLL&McCrajj (Il993h (calcium) 



and Li & McCravl (Il995h (helium). The thermal evolution o f 
the envelope was mod elled in detail by Ide Kool et al ] (119981) . 
and Kozm a & Fransson (1998aib|, in the following KF98a,b) 
computed the spectra from detailed explosion models to study 
the evolution of tempe rature, ionization, and l ine fluxes in the 
200 - 800 day range. iKiaer et al.l (1201 Ol) and iJerkstrand eTall 
(1201 iL J 11 hereafter) analyzed the spectrum in the ^"*Ti -powered 
phase (f > 5 years), including the effects of multi-line radiative 
transfer. 

Modeling of othe r objects includes the work by 
iDessart & Hillieii d2011h . who compared the emergent spectra of 
Type IIP explosion models to SN 1999em, assuming complete 
thermalization of the gamma-rays, but performing a detailed 
radiative transfer calculation. 

In this paper we undertake a detailed analysis of one of the 
brightest and best-observed Type IIP SNe to date, SN 2004et. 
This SN has been subject to progenitor analysis, hydrodynamica l 
modeling, and some qualitative spectral analvsis. lLi et alJd2005l) 
identified a proge nitor candidate in grou nd-based pre-explosion 
images. However. ICrockett et alj d201 ll) showed that this candi- 
date was still visible in 2007 (3 years after explosion), and that 
the source in the original images was a composite of two or three 
sources. They identified an excess flux in the pre-explosion im- 
ages contributed by the true SN progenitor, and for two different 
assumptions about the late-time SN flux determined Mzams - 8 
and 10 M0, with an error of +5/- 1 Mq. With new post-explosion 
imaging and updated methods to determine bolometric correc- 
tions and synthetic colours from model spectra. Eraser et al. 
2012 (in prep.) have revised this to Mzams = H^i Mq. The 
best estimates for luminosity, temperature, and radius are log 
(L/Lo) = 4.8, T = 3600 K, and Rq = 650 Rq. 

Through radiati on-hydrodynamical modeling, 

lUtrobin & Chugail d2009l) derived a large ejecta mass of 
Mej = 24.5 + 1 Mq, the highest value for any Type IIP modeled 
in detail so far They also determined a progenitor radius 
Ro = 1500 + 140 Rq, explosion energy £ = 2.3 + 0.3 B (1 B = 
10^1 ergs) and ^''Ni mass M(5^Ni) = 0.068 + O.OO9M0. 

Spectroscopic analysis is so far limited to some qualitative 
results from the nebular spectrum b v . Sahu et al.. (2006, S06 here- 
after) and MIO. S06 found the [Oi] ^^6300, 6364 luminosity to 
be of comparable to the one in SN 1987 A, and from the similar 
^^Ni mass they conclude that SN 2004et should have a simi- 
lar oxygen mass. Erom this follows a similar p rogenitor mass, 
which was 16 - 22 Mq for SN 1987A dArnett et al. 1989) . MIO 
found the [O i] AA6300, 6364 flux normalized to the ^^Ni mass to 
be 15% lower than in SN 1987 A at 285 days, and assuming all 
other things equal the oxygen mass would then be smaller by the 
same factor With the str ong dependency of oxyge n production 
on progenitor mass (e.g. IWooslev & Weavenll995l) . SN 2004et 
would then be expected to have a progenitor mass only slightly 
lower than SN 1987 A. 

Taken together with the progenitor analysis and light curve 
modeling, there are therefore conflicting results for how massive 
the progenitor of SN 2004et was. The aim in this paper is to 
shed more light on this issue by modeling the photometric and 
spectroscopic nebular phase data. To this end, we use the spectral 
synthesis code described in Jll (see also Maurer e t al. (20li|) 
for testing of the code), with some modifications described in 
Appendix |B] We compare the model spectra produced by using 
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three different explosion models as input (from Mzams =12, 
15 and 19 Mq progenitors), and investigate which produces best 
agreement with observations. This paper is complemented by a 
more g eneral discussion of th e nebular phase spectra of Type IIP 
SNe in lMaguire et al.1 (1201 2l in the following Ml 2). 

2. Observational summary 



SN 2004et was discovered on Sep 27 2004 jZwitter et al.1 



with the explosion date determined to Sep 22 JLi et al 



2004) 



2005) 



The host galaxy was the nearby (5.5 Mpc) starburst galaxy NGC 
6946, which has hosted nine SNe since 1917. The SN was 
classified as a Type IIP, showing broad P -Cygni B aimer lines 
dZwitter et alJ2004l:lF"ilippenko et al.l2004 . Intensive photomet- 
ric and spectroscop ic monitoring was initiat ed, and i s described 
in S06. iMisra etal. (2007. M07 hereafter). iKotak e t al. (2009, 
K09 hereafter) and MIO. The luminosity was among the highest 
of all Type IIP events, both on the 110 day plateau and in the 
radioactive tail phase (MIO). 

T he SN was d etected at radio wavelengths on October 5 
2004 dStockdale et al. 2004.). implying the presence of signifi- 
cant amounts of circumstel lar material. Follow- up observations 
were made with MERLIN ("Bes wick et aDl2004 . GMRT (M07) 
and VLBI ( Marti- Vidal et al. 20071). X-r ay and radio emission 
dRho et alJl2007l: iMartf- Vidal eraiTl2007h indicate the presence 
of a progenitor wind. A mass loss rate of M = 10"^ M© yr ' 
was estimated from the radio turn-on time, assuming a wind 
velocity of v „ - 10 kms~' and a circumstell ar temperature of 
Jcs = 10^ K dChevalier et al . 2006fl whereas iRho et alJ d2007h 
found M - 2x 10"^ M© yr"' from modeling of the X-ray emis- 
sion. 

K09 obtained Spitzer mid-infrared (MIR) spectra at 5 - 35 
fj.m for seven epochs between 295 and 1385 days, as well as 
optical Keck spectra at days 823, 933 and 1 146. The late optical 
spectra appear to be dominated by circumstellar interaction in 
the form of wide (HWZI = 8500 km s'), box-shaped Ha, [Can] 
.17300, and [Oi] ^^6300, 6364, as well as narrow (70kms-i) Ha 
and [Nil] /(/(6548, 6583. Spectroscopic identification of silicate 
dust is reported, with a mass of a few times 10'"^ Mp. Fu rther 
MIR data and analysis is presented by iFabbri etal.ld201l[ Fll 
hereafter). 



2.1. Distance and extinction 

We adopt a distance of 5.5 Mpc f or NGC 6946 , as determined us- 
ing the H I TuUy-Fisher relation (iPiercdl 19941) . and in agreement 
with the a verag e value of several other distance determinations 



(Schmi dt 119931: JSchoniger & Sofud Il994t iKarachentsev et al.l 
[200a:.Sahuet al.ll2006lK 



We adopt an extinction of ^b-v = 0.41, base d on the sum 
of loc al and foreground Nai D absorption lines (IZwitter et al.l 
|2004|) . The g alactic foreground ex tinction towards NGC 6946 is 
Eb-v = 0.34 dSchlegeletal.lll998h . 



2.2. Metallicity 

Metallicity measurements would traditionally come from the 
[O n] and [O m] line flu xes in nebular or H n regions at, or clos e 
to the SN position (e.g.. lSmartt et al.ll2009l; [Modiaz et al.ll201 ll) . 
From the r adial position of SN 2 004et (r^ - 0.92r2 5), and the 
gradient of iPilvugin et al.1 (|2004 . ISmartt et all d2009l) estimated 
the oxygen abundance of the progenitor to 12 H- logO/H = 



8.3 dex, which is a factor of about two lower than solar metal- 
licity. However the determination of absolute oxygen abun- 
dances from the strong line ratio of [On] and [Oiii] has been 
pla gued by large difl ' erence s between the calibrations employed. 
As Moustakas et al.l d2oToh show, the oxygen abundance at the 



The inferred mass-loss rate scales as M oc v„T^ 
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position of SN 2004et determined by two physically motivated 
calibrations is 8.9 and 8.3 dex. This large difference of a fac- 
tor four makes it difficult for us to choose a specific stellar 
metallicity for the pr ogenitor model. The solar value of 8.7 dex 
dAsplund et alJl2009l) is intermediate between the two and as we 
have no clear way of distinguishing between the discrepant re- 
sults we employ a solar metallicity model. We discuss the in- 
fluence of metallicity on the modeling in the discussion (Sect. 
0. 



2.3. Pliotometry 

We use UBVRI photometry presented by S06, M07, MIO and 
Fl 1. S06 obtained photometry for 8 - 541 days, M07 for 14 - 90 
and 370 - 470 days, and MIO for 7-1214 days. Fl 1 obtained 
VRI photometry for 79 - 1412 days. 

For nebular-phase NIR photometry, we use MIO measure- 
ments of J and H magnitudes on days 134 and 141, and JHK 
magnitudes from Fll for days 268, 317, 390, and 652. We also 
use day 307 JHK observations by MIO, which due to their 
importance for the only NIR spectrum taken were recalibrated 
with an extended set of reference stars. The raw data were re- 
duced following standard procedures with the addition of cor- 
rections for quadrant crosstalk and field distortion. We deter- 
mined colour constants from standard fields and approximate S- 
corrections from the standard star spectrum and the filter trans- 
mission curves, to as accurately as possible transform from the 
instrumental system to the standard (2MASS) system. We find 
J = 15.70 + 0.03, H = 15.48 + 0.02, and K = 15.12 + 0.03. 

In the MIR, we use 3.6 //m, 4.5 /im, and 24 /im Spitzer pho- 
tometry from K09. 



2.4. Spectroscopy 

We use a set of sixteen optical, one near-infrared, and three mid- 
infrared spectra from S06, K09, and MIO, a summary of which 
can be seen in Table \T\ We recalibrated the NIR TNG/Amici 
spectrum (day 307) using the new photometry reductions de- 
scribed above. The prism used in these observations has very low 
resolution preventing wavelength calibration using lamp or sky 
lines because of line blending. Instead the observatory provides 
a pixel to wavelength mapping to be calibrated using known 
features in the spectrum by applying a constant pixel shift. We 
used this method to calibrate the spectrum by alignment of the 
strongest lines to the model spectra. The raw data were reduced 
following standard procedures with the addition of corrections 
for quadrant cross talk and the spectrum then extracted and cor- 
rected for telluric absorption using the standard star spectrum. 
The absolute flux calibration was performed by a least squares 
fit of synthetic and measured J and H magnitudes (spectral re- 
sponse in the /T-band region involves significantly higher uncer- 
tainties). 

All spectra displayed in the paper have been cor- 
rected for the host gala xy re cession velocity of 48 
kms ' dSandage & Ta mmannI Il98lh . The optical and NIR 
spectra have been c orrected for ex t inction with ^b-v - 0.41 
(Sect. IZTT ) and the ICardeUi et aP dl989h extinction law. No 
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correction was applied to the MIR data, for which the extinction 
is negligible. 



Table 1. Summary of spectroscopic observations used in this pa- 
per. 



Phase 


Date 


Wl. coverage 


Res. 


Ref. 


163 d 


2005-03-03 




- 0900 A 


7 A 


S06 


169 


2005-03-09 


JO\J\J 


- 7730 A 


24 A 


MIO 


212 


2005-04-21 




- o?fln A 


7 A 


S06 


249 


2005-05-28 


J~J\J\J 


- 0900 A 


7 A 


S06 


257 


2005-06-05 


JOUU 




13 A 


MIO 


284 


2005-07-02 


^uuu 




36 A 


MIO 


295 


2005-07-13 


5.2 ■ 


— 31 /jm 


~ 0. 1 ■ /l^,^ ^m 


K09 


301 


2005-07-19 


3500 


- 9200 A 


7 A 


S06 


307 


2005-07-25 


0.8- 


2.54 yum 


0.2 • /l^m ^im 


MIO 


314 


2005-08-01 


3500 


- 9200 A 


7 A 


S06 


341 


2005-08-28 


3400 


- 8000 A 


13 A 


MIO 


350 


2005-09-06 


5.2- 


- 35 yum 


~ 0.1 • /l^m//m 
24 A 


K09 


384 


2005-10-10 


3500 


- 7730 A 


MIO 


391 


2005-10-17 


3500 


- 9200 A 


7 A 


S06 


401 


2005-10-27 


3800 


- 7730 A 


25 A 


MIO 


401 


2005-10-27 


3500 


- 9200 A 


7 A 


S06 


408 


2005-11-03 


4000 


- 7730 A 


25 A 


MIO 


428 


2005-11-23 


3500 


- 9200 A 


7 A 


S06 


451 


2005-12-16 


5.2- 


- 20 yum 


~ 0.1 • /l^„,//m 


K09 


465 


2005-12-30 


3500 


- 9200 A 


7 A 


S06 



3. Modeling 

The spectral synthesis code we use is described in detail in 
Jll. In summary, we compute the energy deposition by ra- 
dioactive isotopes by performing an effective opacity trans- 
fer of the gamma-rays, and assuming on-the-spot deposition 
for the leptons and X-rays. The relative fractions of the en- 
ergy going into heating, ionization, and excitation are cal- 
culated as described in iKozma & FranssonI (Il992h . and the 
equations of statistical and thermal equilibrium (including 
non-thermal heating/ionization/excitation, collisional (thermal) 
excitation/de-excitation, radiative excitation/de-excitation, pho- 
toionization, radiative recombination, charge transfer, photo- 
electric heating, free-free heating, line cooling, recombination 
cooling, and free-free cooling) are then iterated with a Monte 
Carlo simulation of the radiation field until convergence has 
been achieved. The radiative transfer treatment uses the Sobolev 
approximation. The method is a quasi-A-iteration, but differs 
in that it treats continuum and part of the line scattering on- 
the-fly in each radiation field simulation, rather than including 
them in the emissivity functions. This improves the convergence 
properties by allowing photons to travel longer distances than 
one mean-free-path in each iteration. The emergent spectrum 
is obtained by binning the escaping Monte Carlo packets, and 
smoothing with a Gaussian of FWHM = 600 km s"^ . 

The code has been updated in some respects compared to 
the version used in Jll. Because the computations here are for 
earlier epochs than in Jl 1, photoexcitation/de-excitation rates in 
the NLTE solutions cannot be approximated to zero. We there- 
fore implemented a modified Monte Carlo operator to allow for 
partial deposition of photo-excitation energy along the photon 
trajectories (we keep the coupling partial to avoid approaching a 
classic A-iteration, and keep the NLTE solutions limited in size). 



The details of the new treatment can be found in Appendix iBl 
We have also added electron scattering, as well as the influence 
of line overlap on Lya and Ly^S Sobolev escape probabilities, 
also detailed in Appendix iBl 

3. 1 . Explosion model and zoning 

As input we use t he 12 Mp, 15 Mg , and 19 M© explosion 
models computed bv lWooslev & He^ed (l2007l WH07 hereafter), 
which are non-rotating stars of solar metallicity evolved with the 
KEPLER hydrodynamics code and then exploded with a piston, 
endowing a final kinetic energy of 1.2 B to the ejecta. The mass- 
loss prescription in these models is based on the empirically es- 
tablished relations by Nieuwenhuijzen & de Jager ( 1990), both 
for the main sequence and the RSG phase. The mass loss for the 
three stars are 1.2, 2.3, and 3.8 M©. 

We divide the ejecta into eight distinct zone types : Fe/He, 
Si/S, O/Si/S, O/Ne/Mg, O/C, He/C, He/N, and H, named after 
their dominant constituents. The mass of the Fe/He zone is set to 
give a total ^^Ni-mass of 0.062 M© in the ejecta (see Sec. 13.3b . 
Table |2] shows the zone masses in the three models we use, and 
Tables lAniA.3l in the appendix show their detailed composition. 



Table 2. Zone masses (in Mq) of the three different explosion 
models used. 



Zone / Model 


12 


15 


19 


Fe/He 


0.10 


0.087 


0.076 


Si/S 


0.061 


0.080 


0.16 


0/Si/S 


0.13 


0.24 


0.17 


O/Ne/Mg 


0.14 


0.45 


1.9 


O/C 


0.16 


0.43 


0.60 


He/C 


0.16 


0.74 


0.88 


He/N 


0.92 


0.43 


0.40 


H 


7.7 


8.4 


9.5 


Total 


9.3 


10.9 


13.6 



The explosion models are ID, and therefore lack the macro- 
scopi c mixing which is kn own to occur (e.g., iKifonidis et al.l 
120061: iHammer et al.llIoToh . We therefore distribute the zones 
over a core, where each zone is distributed as A^ci spherical 
clumps, completely macroscopically mixed in velocity, and an 
envelope, where the zones are concentric shells, with no mix- 
ing. The Fe/He, Si/S, O/Si/S, O/Ne/Mg, O/C, and He/C zones 
exist only in the core, whereas the He/C, He/N, and H zones ex- 
ist partially in the core (fractions xue for the He zones and xh 
for the H zone), and partially in the envelope (fractions 1 - jchb 
and 1 - jch). The first two shells in the envelope are the unmixed 
He/C and He/N zones. External to these zones are logarithmi- 
cally spaced {V-JVi-i - 1.2) shells of the H-zone material. The 
velocities of the two unmixed He-shells are set so their densities 
equal that of the first hydrogen shell. 

The density of the H-envelope is set to follow the density 
profile of the input model, but scaled down with a constant to 
account for the mixing of material into into the core. This down- 
scaling was about 10% for all three models. The envelope is ter- 
minated at lO'* kms"', as the density decreases steeply beyond 
a few thousand kms and little line or continuum opacity re- 
mains at higher velocities than this in the nebular phase. 

Although our artificial mixing alters the velocity field for 
each species compared to the ID input models, the explosion 
energy of 1.2 B is preserved to within 3% in all models. Fig.[T] 
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shows the velocity field for the 15 Mq model that we use, com- 
pared to the ID hydro simulation. 
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Fig. 1. The density profile in the ID hydrodynamical model 
(solid lines) and in the artificially mixed model used in this pa- 
per (dashed lines). Note that in the artificially mixed model, the 
zones in the core are distributed as clumps with individual filling 
factors, the densities plotted are the densities within each clump. 



3.2. Model parameters 

The parameters of the model are the velocity dividing the core 
and the envelope (Vcore), the fractions of the He and H zones 
mixed into the core (xHe and xh), the filling factors of the differ- 
ent zones in the core (/pe/He, /si/s, /o/si/s, /o/Ne/Mg, fo/c, /hb/c, 
fiie/N, /h), the number of clumps per zone in the core (A^ci), and 
the dust optical depth function (Td(f)), which we assume to be 
gray. Most of these parameters can be constrained from observa- 
tions or theoretical arguments, as described below. 



3.2.1. Core velocity Vcore 

The core velocity is to first order found by inspecting the line 
widths of [Oi] /1/16300, 6364 and [Feii] ^7155, which were 
mainly emitted from the core in SN 1987A (KF98b). The wings 
of these lines may, however, be produced by primordial oxy- 
gen and iron in the envelope (Ml 2), and are therefore unsuit- 
able markers of the core width. A more relevant starting point 
is the half-width-at-half-maximum, HWHM, which is ~1200 
kms-i for the [Oi] AA6300, 6364 lines (M12). Inspection of 
the [Fell] /17155 line gives similar values. The emission line 
profile from a homogeneous sphere expanding homologously is 
parabolic with Vhwhm = 0.71 Vcore, so Vcore can be estimated 
to 1700 kms"'. We find that a value Vcore - 1800 kms"' best 
reproduces the line profiles, and is chosen for the modeling. 

3.2.2. In-mixing fractions xue and xh. 

Several 2D and 3D hydrodynamical calculations have shown 
that 1-2 Mq of the H-envelope gets mixed with core mate- 
rial in the explosion (HeraiTt&J3enz 199 1 ; Mueller et al. 199 f[ 
Herant & WooslevI 119941: iKifonidls et al.l 2006; H ammer et al.l 
20101) . This was also found from modeling of the Ha evolu- 



tion in SN 1987A (KF98b'). lHerant & WooslevI (1 1994 obtain 1.4 
Mo of H (20% of the total) and 1.3 Mg of He (85% of the to- 
tal) inside 1800 kms for the 2D explosion of a 12 Mq RSG. 
The 3D explosion sim ulation of a 15.5 M© blue supergiant by 
iHammeret all (1201 Ol) shows 60% of the He and and 15% of the 
H mixed inside Vcore = 1800 kms Based on this we adopt 
xue - 0.60 and xh = 0.15. 

3.2.3. Filling factors 

The filling factor of the iron zone is determined by the energy 
input from the ^''Ni deca y during the first weeks, creating an 
"iron bubble " in the core (IWooslevl[T98l iHerant & Benzll 19911: 
lBaskolll994 . As we show in Sect. 14.41 analysis of the Fe/Ni/Co 
fine-structure lines can be used to determine the filling factor of 
this zone. We find /pe/He ~ 0.15, which is used throughout the 
rest of the paper 

The oxygen density can be derived from the evolution of the 
[O i] /16300/6364 line ratio, which is found to be similar to the 
one in SN 1987A (M12). We therefore use the same value for the 
oxygen de nsity as derived for SN 1987A, rip = 6.2x10'° c m~^ at 
100 days dSpvromilio & Pintolll99U iLi & McCravlll992h . This 
determines /o/si/s, /o/Ne/Mg, and /o/c, giving a scahng propor- 
tional to the oxygen mass. 

The Si/S zone was found by Jl 1 to have expanded to a low 
density in SN 1987 A, and is given a density ten times lower than 
the oxygen zone density. 

The remaining core volume is divided between the mixed- 
in He/C, He/N, and H components to give them equal number 
densities. 

The resulting filling factors for the three models are listed in 
Tabled 



Table 3. Filling factors used for the core zones the three models. 



Filling factor \ Model 


12 Mo 




15 Mo 




19 Mo 




,/fc/Hc 


0.15 




0.15 




0.15 




,/si/s 


1.6 X 10 


-2 


2.2 X 10- 


2 


4.4 X 10- 


-2 


/o/si/s 


7.9 X 10 


-3 


1.5 X 10- 


2 


5.8 X 10- 


-3 


/o/Ne/Mg 


6.8 X 10 


-3 


2.3 X 10- 


2 


0.10 




/o/c 


8.7 X 10- 


-3 


2.4 X 10- 


2 


2.8 X 10- 


-2 


/hc/C 


1.4 X 10- 


-2 


7.3 X 10- 


2 


7.2 X 10- 


-2 


/hc/N 


0.10 




4.4 X 10- 


2 


3.4 X 10- 


-2 


/h 


0.69 




0.65 




0.57 





3.2.4. Number of clumps A^ci 

Several arguments suggest that the number of clum ps pro- 
duced in SN explosions is high. For SN 1987 A, IChugail 
d 19941) could reproduce the statistical fluctuations in the [Oi] 
AA63QQ, 6364 lines with 2000 clumps. Multi-dimensional mod- 
els i ndicate strong mixing and a high degree of fragmenta- 
tion dKifonidis et al.ll2006t iHammer et al.ll2010h . The strong Oi 
/11.129 ;um line observed in SN 1987A requires synthesizecQ 
oxygen to be mixed with hydrogen down to very small scales 
in order for the Bowen fluo rescence to occur before Ly/3 is con- 
verted to Paa (l01ivalll9^3l) . also supporting a large number of 



^ We use the word "synthesized" for material made during either hy- 
drostatic or explosive burning. Oxygen is made mainly by hydrostatic 
burning. 
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clumps. Based on these results, we choose the number of clumps 
A^ci = 10^ for our modeling here. 

3.2.5. Dust 

Both S06 and MIO note that the [O i] /16300 and Ha lines show 
shifts to the blue at around 300 days, and conclude that dust 
formed at this time. We have extended this analysis of the dust 
formation process by measuring the position of the line peaks 
forNai/158943, [Oi] ^6300, Ha /I6563, [Fen] ^7155 and [Can] 
/I730fl for all epochs in the 163-465 day range. 

We measured the positions of the peaks in two ways; first 
by determining the wavelength of the maximum (within ± 2000 
km s"' of the rest wavelength), and secondly by fitting a second- 
order polynomial to the part of the line that is within 15% of 
the peak flux. The two methods gave similar results. The result- 
ing measurements with the second method are plotted in Fig.|2] 
together with straight-line fits. 
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Fig. 2. Velocity shifts of the peaks of the strongest emission lines 
in the observed spectra as function of time. The symbols are es- 
timates of the peak positions from fitting a parabola to the part of 
the line with a flux higher than 85% of the peak value (95% for 
[Fen] /17155). The lines are linear fits to the measured values. 



All emission lines show a continuous blue-shifting trend 
with time. However, all lines, except [Oi] /16300, show a red- 
shift of a few hundred kms ' at early epochs. One possibility 
is that the SN has a higher recession velocity than the galactic 
one of 45 kms However, positioned far out from the cen- 
ter in a galaxy seen almost face-on, a peculiar velocity of sev- 
eral hundred kilometers per second is unlikely. Another pos- 
sibility is that these redshifts are due to electron scattering in 
the ejecta, wh ose opti cal depth is of or der unity i n the early 
nebular phase (IChugaCl977MDessart & Hilliej|201 Ih . However, 



We find the Nai /1/15890, 5896 emission to be much stronger than 
He I /i5876 in all our models. The Na i lines are optically thick through- 
out most of the ejecta for the epochs studied here, so Nai /i5890 emis- 
sion will scatter in Nai /i5896, making 5896 A the suitable wavelength 
choice. 

' This is the peak wavelength if the lines are optically thin (as ex- 
pected for the primordial calcium dominating the emission (M12)) and 
have Gaussian shapes with Vfwhm = 1300 kms"'. 



while electron scat tering will produce a dis t inct red wing 
in the line profile (Auer & van Blerkom 1972'; 'Chugaf '1977t 
Ifransson & Chevalier 1989), Fransson & Chevalier (1989) did 
not find any redshift of the peak in their simulations. We there- 
fore favor a third possibility, that a fundamental asymmetry of 
the ejecta is the explanation. 

Whatever the reason is for the early red-shifts, later on the 
line peaks change into blue-shifts, which is most likely caused 
by dust. The average value of the slope of the linear fits is -2.0 
kms^V day. Denoting the velocity shift as function of time by 
AV(t), and taking the dust formation to begin at day 250, we 
then write 



AV(t) = -2.0 X (f - 250 d) km s"', t> 250 d . 



(1) 



The dust optical depth Td(f) can then be estimated from the 
formul a for a uniformly emitting and absorbing sphere of veloc- 
ity Vi dLucv et al.lll991h : 



Am 

Vd 



= -1 + 



ln[l+Td(0] 

Td(f) 



(2) 



The solution for T^it), for Vd - 1800 km s ' , is plotted in Fig. |3] 
(red dashed line). 



0.8 
0.7 
0.6 
0.5 
H°0.4 
0.3 
0.2h 
0.1 

ft 



Fit to line peak blushifts 
Fit to MIR dust emission 
Function used for the modeling 



50 



200 



250 



300 350 
Time [days] 



400 



450 



500 



Fig. 3. The (gray) dust optical depth function Td(f) derived from 
the observed dust luminosity (blue, solid), from the blue-shifts of 
line profiles (red, dashed), and the one we use for our modeling 
here (black, dot-dashed). 

The dust hypothesis is supported by Spitzer observations of 
a warm MIR-continuum at 500 K (K09, Fll). In K09, a triple- 
component black-body fit reveals a steadily increasing fraction 
of the ^^Co-deposition being re-radiated by a warm component, 
which is attributed to dust formed in the ejecta. By equating 



a(t) 



Ldep(t) 



= 1 - pit) 



(3) 



where Lwaim(0 is the luminosity of the warm component, 
Ldepit) is the energy deposition, and p{t) is the escape proba- 
bility from a uniformly emi tting sphere with constant opacity 
dOsterbrock & Ferlandll2006l) . 



Pit) 



4rd(f) 



1 



1 



1 



1 



2Td(f)2 \Tdit) 2rd(f)- 



,-2tj(0 



(4) 



6 



Anders Jerkstrand et al.: The progenitor mass of the Type IIP supernova SN 2004et from late-time spectral modeling 



the dust optical depth function Td(0 can be determined. The re- 
suhing values (with L„arm(0 and Ldep(0 taken from K09) are 
plotted in Fig. [3] (blue, solid line)0, and are compared with the 
values derived from the line-shifts. Within the errors expected 
from the uniform sphere assumption, the independent deriva- 
tions agree reasonably well. 

For our modeling, we employ a functional form for the dust 
optical depth lying between the Spitzer and line-shift fits, given 
by 

Td(f) = 1-8 X 10"^(f-250d), f> 250 days, (5) 

shown as the black dot-dashed line in Fig. [3] The extrapolation 
after 500 days is uncertain, but is not relevant for any of the 
conclusions in this paper. 

3.3. ^^Ni mass 

The ^^Ni mass can be determined by comparing the estimated 
bolometric luminosity in the early tail-phase with the theoretical 
value of fully trapped ""^Co deposition, which is given by 

L3.co(0 - 9.92 X lO-*' (e-'"''-' ' - e-"''-' erg s"' .(6) 

^° 0.07 Ma ^ ; 5 V / 

Apart from the assumption of full trapping, one also has to as- 
sume that the deposited energy is instantaneously re-emitted, 
and that no other energy source has any influence. 

We used the UBVRIH photometry from day 142 (MIO), 
together with an extrapolation outside the //-band (see below) 
to estimate the bolometric luminosity at this epoch. The lin- 
early interpolated, dereddened UBVRIH flux at 142 days is 
1.82 X 10"*' ergs '. The UB VRIJHKLMN^NjN^Qn photome- 
try of SN 1987 A at day 135 dHamuv et al.ll 19881: iBouchet et all 
119891) shows 26% of the total (linearly interpolated) U - Qo flux 
to be outside the U - H range. Assuming a fraction of 0.26 
also for SN 2004et, the extrapolated bolometric U - 2() fux is 
2.46 X 10"*' erg s"' . From Eq. (|6]l, this corresponds to a ''''Ni mass 
of 0.062 Mq. This is in good agreement with values derived by 
S06, M07, and MIO. 

3.4. Molecules 

There are now a handful of detections of CO in core-collapse 
SNe, s ugges ting that its formation is common (e.g. lGerardv et al.l 
120001 l2002h . In SN 1987A, CO emission was detected al- 
ready in the p hotospheric p hase, and persisted to o ver 600 days 
dSpvromilio et al.„1988-.Bouchet & Danzigeiill993h . 

Also SiO was formed in SN 1987A, detected fro i n day 
~250 an d persis ting to over 600 days dRoche et al.l Il99lt 
iBouchet & Danzi gerlll993h . It has fewer detections in other SNe 
compared to CO, likely because it lacks a strong spectral feature 
in the NIR. 

The 2.3 fim first overtone band of CO was detected in SN 
2004et in the NIR spectrum taken at 307 days (MIO). No earlier 
nebular NIR spectra exist to determine exactly when this feature 
emerged. 

The 8 fim fundamental band of SiO was detected in the first 
Spitzer spectrum at day 295 (K09). It could then be followed un- 
til at least day 450, after which it became difficult to distinguish. 
No spectra covering the fundamental band of CO (4.6 //m) were 

* We only include the results up to day 465, since after that time 
dust emission from circumstellar interaction affected the MIR emission 
(K09). 



taken, but its presence is revealed by the photometry (see Sect. 

ESS. 

Models show that CO forms in the O/C zone, and SiO 
forms in the O/ Si/S zone, with little fo rmation of either in the 
O/Ne/Mg zone (lLiu&Dalgarndll995h . Our models do not in- 
clude molecules, and because these l ikely determine the tem per- 
ature of the O/C and O/Si/S zones dLiu & Dalg arno 1995*). we 
need some prescription for how to set these temperatures. 

ISpvromilio et al ] dl988l) found, assuming LTE and optically 
thin conditions, the temperature of the CO clumps in SN 1987 A 
to decfine from 3000 K to 1200 K between 192 and 350 days. 
iLiu et al.l d 19921) found that inclusion of optical depth and NLTE 
effects gave a temperature evolution of 4000 - 1900 K. We as- 
sume that the temperature in the CO zone follows a linear inter- 
polation of this result, but levels off' at 1000 K at late times 

ro/c-zone(0 = max{4700 - 13 • (r - 140 d) K, 1000 K) . (7) 

For SiO, iLiu & Dalgarnd d 19941) found a temperature evolution 
from 2000 K to 1000 K between 257 and 519 days in SN 1987 A. 
We use the same interpolation method here 

ro/si/s-zo„e(f) = max{2400 - 3.8 • (f - 140 d) K, 1000 K) . (8) 

Our results are not sensitive to the exact prescription for the tem- 
perature of these zones as they have too small mass and are too 
cold (assuming that some level of molecular cooling occurs) for 
any significant atomic thermal emission. 



4. Results 

4.1. Temperature evolution 

For later guidance in the interpretation of line ffuxes. Fig. |4] 
shows the model temperature evolution for the diff'erent zones, 
for the 15 Mq model. These results may also be of use in in- 
terpreting thermal emission line luminosities in other SNe with 
similar ^^Ni mass. 
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Fig. 4. Evolution of model temperatures in the various core 
zones in the 15 M© model. 



7 



Anders Jerkstrand et al.: The progenitor mass of the Type IIP supernova SN 2004et from late-time spectral modeling 




13 






7 14 




— ^ — ^ — \ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — 


1 12 


^ — 1 — 1 — 1 — 1 — 1 




1 13 

-D 




— ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — \ — ^ — ^ — ^ — ^ — \ — ^ — I|"'T - 


T3 
CM 

o 12 


1 1 1 1 1 




' 13 

X 






CM 


^ \ 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


7 13 







300 400 500 600 700 
Time [days] 



100 200 300 400 500 
Time [days] 



600 



700 



Fig. 5. Normalized photometric evolution (all magnitudes have been subtracted by f/102 d to normalize to the ^''Co decay) of the 
(reddened) 12 (blue, dot-dashed), 15 Mq (black, solid), and 19 (green, dashed) models, compared to observations (red 
points). T.P. means 'this paper'. The MlOs, S06s, and TPs data points are synthetic photometry from the observed spectra. 



4.2. Photometric evolution 

Fig. |5] shows the observed (normalized) U,B,V,R,I,J,H and 
K photometric evolution compared to the models. For clarity 
we have divided all fluxes with the exponential decay factor 
for ^^Co, which in magnitudes corresponds to a subtraction by 
f/98 d. This reduces the span of magnitudes from ~8 to ~2, al- 
lowing greater detail to be seen. 

In addition to the observed photometry, we also plot photom- 
etry obtained by convolving the observed spectra with the stan- 
dard filter functions (filled (black) symbols). As is clear from 
the 7? -band and /-band plots, significant differences may occur 
here, due to the difficulty of transforming from telescope filters 
to standard filters for sources with strong emission lines, such as 
nebular phase SNe. 

The models diff'er little in the [/-band, all agreeing well with 
the data. The agreement (to within 20%) is quite striking, con- 
sidering the enormous complexity of the radiative transfer in 
the UV. The relative brightening with time (up to the on-set 
of gamma-ray leakage at ~400 days) is due to decreasing line- 
blocking. 

In the B-band, the models are also similar, reproducing the 
general time-evolution well. As in the U -band, the relative flux 
initially increases as line blocking decreases in efficiency. The 
small differences between the models in both U and B can be 
understood as a consequence of most of the radiation here origi- 
nating from far out in the envelope (which has a similar structure 
in all models), the core being highly opaque due to line blocking. 



From the V-band and on, line blocking is less efficient, and 
differences due to the different cores start to emerge. In V, the 
12 and 15 Mq models are in good agreement with the observa- 
tions at all times, whereas the 19 model is too bright. The 
y-band includes Na i AA5890, 5896 and [O i] AA6300, 6364, and 
it is mainly too strong emission in these that cause the overpro- 
duction in the high-mass model (Sect. l43] l. 

The /?-band behaves similarly to the V-band. The band is 
dominated by Ho-, but is also influenced by [O i] AA630Q, 6364, 
especially at late times. The fluxes in both lines are too strong in 
the 19 Mo model (see below), overproducing the /?-band. The 12 
and 15 Mq models are better, but still show some discrepancies 
as they also over-produce Ha (Sect. l431 l. 

The /-band shows significant differences between the ob- 
served photometry and the synthetic photometry from the ob- 
served spectra, up to half a magnitude during the first year. The 
/-band is influenced by the [Can] 7291, 7323 doublet as well as 
the Ca II IR triplet, both lying on the edge of the filter As men- 
tioned before, this implies that SN /-band photometry has to be 
used with considerable caution. 

In the NIR, the / and H data points at 134 and 141 days are 
important in showing that the models reproduce the flux reason- 
ably well also in the NIR in the full gamma-ray trapping phase, 
which together with the good agreement in the other bands pro- 
vide confirmation that the ^*Ni mass is well determined. For the 
/iT-band, one should be aware that our models do not include the 
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Fig. 6. Optical (dereddened) spectrum at 401 days (red) compared to the 12 (blue) and 19 M© (green) models. 
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Fig. 7. Observed (dereddened) spectra (red) and 15 Mq model spectra at 212, 401, 
smoothed with a Gaussian of FWHM = 600 km s 'to improve the S/N. 



and 465 days (black). All spectra have been 



CO first overtone at 2.3 fim (which lies at the edge of the band), 
so we expect the models to be on the low side. 

It is clear that the progenitor mass has its strongest influence 
on the V, R, and H bands, which are therefore valuable observ- 
ables in the nebular phase. In U, B, I, J, and K, the differences 
are smaller. 



4.3. Spectral evolution 
4.3.1. Optical 

Fig. |6] shows the optical model spectra for the 12 and 19 
models compared to the observed (dereddened) spectrum at 401 
days. It is clear that the main influence of the progenitor mass 
are the Na i ^^5890, 5896 and [O i] AA63Q0, 6364 lines. The ob- 
served spectrum shows fluxes in these lines that fall in between 
the 12 Mq and 19 Mq models. 



Also Ha, [Can] AA7291, 7323, and the Can IR triplet fluxes 
differ between the models, although by less. These lines are 
mainly formed in the hydrogen zone, and are therefore not di- 
rect indicators of the core mass, but rather depend on the quite 
uncertain pre-SN mass loss and the explosive mixing. 

The spectrum below 6000 A is a complex mix of mainly 
iron-group lines formed by scattering and fluorescence in the hy- 
drogen envelope. Because of the similar envelopes in the mod- 
els, the model spectra differ by little in this spectral region. One 
exception is the Mgi] /14571 line, which arises from the metal 
core. 

The 15 Mq model spectra at 212, 401, and 465 days are plot- 
ted in Fig. [T] We see that this model is generally successful in 
reproducing the observed spectra, and in particular reproduces 
the [Oi] AA630Q, 6364 lines well at all epochs. This is further 
corroborated by Fig. |8] where we plot the evolution of Mg i] 
i4571, Nai AA5890, 5896, [Oi] AA63Q0, 6364, [Neii] ^12.81 
fim, Ha, [Fen] ^^7155,7172, [Can] ^^7291,7323, and the Can 
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Fig. 8. Observed (dereddened) line fluxes (red points), normalized to the bolometric flux of 0.062 of ^^Co at a distance of 5.5 
Mpc, compared to model values (12 : black, solid, 15 Mq : blue, dot-dashed, 19 ; green, dashed). See the text (Sect. l43?T] ) 
for details on the line flux measurements. 



IR triplet, from the observed spectra as well as from the three 
different models. Both observed and modeled fluxes have been 
normalized to the bolometric flux of 0.062 of ^^Co, assum- 
ing a distance of 5.5 Mpc. The resulting numbers therefore rep- 
resent the fraction of the radioactive decay energy re-emerging 
in the various lines. The optical line flux measurements were 
done by integrating the spectra over V^ut = +2700 kms"' and 
subtracting the continuum level, which was taken as the linear 
interpolation of the minimum flux values on the right and left- 
hand side, within Vcont - ±l-25ycut = +3375 kms"'. These 
velocity limits were found to give reasonable extractions of the 
line fluxes for all lines and epochs, except for Ha for which we 
used y^ut = ±7200kms-i and Vcont = 1-25^01,1 = ±9000 km s"' 
to include the broad wings. Because we perform the same opera- 
tion on model and observed spectra, the choice is not critical for 
the comparison. 

Diagnostic lines of the core mass The most important lines 
as diagnostics of the core mass are [Oi] /1/16300, 6364. These 
lines are distinct at most epochs, relatively undisturbed by line- 
blending and line blocking (the re is, however, some iro n line 
emission blended with them, see lDessart&Hime3(l201lh ). and 
have at least half their flux coming from synthesized oxygen 
(M12). In addition, oxygen is mainly neutral in the core, while 
magnesium and sodium are mainly ionized, meaning that the 
thermal emission by O i is quite insensitive to the ionization bal- 
ance. In all models, the thermal emission by [Oi] AA63QQ, 6364 
from the 0/Ne/Mg zone corresponds to 50 - 60% of the heat- 
ing of this zone, at all times after 300 days. Combined with the 
constraints on the density that the 6300/6364 flux ratio provides 



(Ml 2), the doublet lines are excellent diagnostics of the core 
mass. The main complication is that primordial oxygen in the 
hydrogen zone contributes to the emission as well (M12). In the 
15 and 19 M© models, however, there is so much synthesized 
oxygen that this contribution is quite small. 

Also the Mg i] AA51 1 line is a potentially useful diagnostic of 
the O/Ne/Mg zone. According to the models, this is the strongest 
optical Mgi line, followed by Mgi ^/15167, 5172, 5183, which 
is too weak to be detected. The Mgi] /14571 line is, however, 
also rather weak and difficult to extract. The first clear detection 
is in the spectrum taken at 284 days. Before that, complex line 
blending and blocking makes it difficult to distinguish, and even 
later it is possible that the emission is blended. The last spectrum 
at day 465 has too poor S/N to assess the line. 

None of the models produce any distinct Mgi] AA51\ un- 
til 200 - 300 days after explosion, in agreement with the ob- 
servations. The evolution after that is in best agreement with 
the 15 Mq model, but the model differences are quite small. 
Observations after 500 days would have allowed a more con- 
clusive comparison between the models to have been made with 
regard to this line. 

Mg II is the dominant ion in the O/Ne/Mg zone at all epochs, 
so the Mgi lines have strong contributions by recombination. 
The transition has a thermal emission of about 10% of the [Oi] 
/1/16300, 6364 emission. This is not enough to account for the lu- 
minosity of the line (Fig. [8]), so recombinations indeed dominate 
the emission. 

The Na i D resonance lines are optically thick quite far out 
in the H-envelope (to 6200 kms ' at day 300). The lines are 
therefore formed by a combination of P-Cygni scattering and by 
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Fig. 9. Near-infrared (dereddened) spectrum at 307 days (red) compared to the 12 Mq (blue), 15 M© (black), and 19 (green) 
models. The boxed regions at 1 .5 - 1 .2 jum and 1 .4 - 1 .7 fim show particular sensitivity to the core mass. Note that CO emission is 
not included in the model spectra. 



emission by synthesized sodium in the core. The D lines are re- 
sponsible for some of the cooling of the O/Ne/Mg zone, but most 
of the emission is recombination. The observed flux is in good 
agreement with the 12 and 15 models, but is strongly over- 
produced by the 19 model. Less influenced by line block- 
ing and blending than the Mg i] /14571 line, but still sensitive to 
charge transfer (Jll) as well as the hydrogen envelope proper- 
ties (through the scattering), this line is intermediate to the Mg i] 
A451 1 and [O i] AA630Q, 6364 lines in model parameter sensitiv- 
ity. 

We do not reproduce the [Neii] 12.81 fim line well with any 
of the models. We discuss this line more in the section about the 
MIR lines. 



Diagnostic lines of tine liydrogen and lielium zones In addi- 
tion to the the hydrogen and helium lines, also many metal lines 
are dominated by emission from the hydrogen and helium zones. 
The strongest lines from these zones are discussed in detail in 
Ml 2, and we will therefore here only make some specific points 
relevant for this paper. 

The hydrogen envelope mass is similar in the three models, 
7.7 - 9.5 Mq (TableO. We also assume the fraction mixed into 
the core to be the same in all models (xn - 0.15). There are 
therefore no major differences in the H-zone lines between the 
models. Still, as Fig. [8] shows, the lines from this zone do dif- 
fer somewhat, as the hydrogen and helium zones have different 
densities for their in-mixed components, and are also exposed 
to different amounts of gamma-rays as well as diffuse radiation, 
depending on the metal content of the core. 



Ha is formed by recombinations following photoionizations 
in the Balmer continuum up to ~500 days, and recombinations 
following non-thermal ionizations from the ground state after 
that. With less core material to absorb the gamma-rays in the 
lower-mass models, the Ha flux is somewhat higher, although 
the hydrogen mass is smaller. The observed flux is lower than 
in all models. As most of the model flux comes from the in- 
mixed hydrogen component, a reduction in in-mixing below 
15% would likely improve the fits. Note, however, that the Ha 
flux in SN 1987 A was about a factor of two weaker than ex- 
pected (from the oth er hydrogen re combination lines) in the 
200 - 500 day range dXu et al.lll992h . still without an explana- 
tion. 

The [Ca ii] A A7291, 7323 lines ar e significant coolants of the 
hydrogen zone dLi & McCravl [19931 KF98b, M12). The mod- 
els differ by little, both agreeing reasonably well with the data. 
The ratio of the Can triplet lines to th e 7291, 7323 A lines 
depends on the temperature and densi ty (Fransson & ChevaUetl 
119891: iFerland & Perssonlll989t iLi & McCrav.,1993i) . Our mod- 
els give a slightly too high triplet to 7291, 7323 A ratio, which 
indicates that our density in the Ca ii emitting gas may be too 
high. 



4.3.2. Near-infrared 

Fig.|9]shows the day 307 NIR spectrum compared to the models. 
It is clear that the core mass has a significant influence on the 
spectral features at 1 . 15 - 1 .2 jum and 1 .4 - 1 .7 fim. The plateau- 
like feature at 1.15 - 1.2 jum is a blend of Ci /II. 176, 1.181 
jum (C/O zone + He/C zone), Mgi .11.183 jum (O/Ne/Mg zone). 
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Si I /II. 200 /im (Si/S zone + O/Si/S zone), as well as several Fe i 
lines from both the Fe/He zone and from the O/Ne/Mg zone. The 
1 .4- 1 .7 yum range contains emission by C i /II .454 fim (C/O zone 
+ He/C zone), Mgi /II. 504 jum (O/Ne/Mg zone), [Sii] /II. 607, 
1.645 fim (Si/S zone), and [Feii] /II. 534 ;um, 1.600 fim, 1.644 
fj.m, 1.681 fim (Fe/He zone). As an important complement to the 
optical analysis, we find that the 15 model produces better 
fits than the 12 and 19 M© models for these metal emission lines 
as well. 

The Oi + Nai feature at 1.13 jum may seem to suggest a 
higher-mass progenitor than 19 Mq. However, the O i /II . 129 /im 
line (which dominates the blend) is produce d by H i LyP b ranch- 
ing into O I /1 1025. 76 2pVP)-3d{^D"), and lOlival ([T993b found 
that much of this branching must occur over zone boundaries, 
with Lyy6 photons scattering into regions of synthesized oxygen 
while trapped between n - I and n = 3 in H i. Our modeling 
does not treat such line overlap between zones (although we do 
treat the overlap within each zone, see App. IB. 51 ), and we can 
therefore not properly model the line. 

The strongest feature in the NIR is Hei /1 1.083 fim + [Si] 
/II. 082 fim, which is a bit too strong in the models. Of the hy- 
drogen lines, both Paff and Pa/3 are also overproduced. Together 
with the overproduction of Ha, this suggests that the He and 
H content of the models may be too large, and/or too strongly 
mixed into the radioactive core. 

4.3.3. Mid-infrared 

As shown by K09, and discussed in Sect. 13.21 the mid-lR is 
dominated by a combination of warm dust emission from the 
ejecta and cold dust emission from the circumstellar environ- 
ment, with superimposed atomic and molecular emission lines 
from the ejecta. To reproduce the spectrum in this region it is 
therefore necessary to include these dust components. For this 
we repeat the modeling in K09, although with slightly different 
parameters, since we here include also the detailed contribution 
by atomic emission (which was approximated as a hot blackbody 
in K09). 

For ejecta dust extinction we use the smoothed Draine & 
Lee model with "as tronomical" silicate dPraine & Led [19841 
iLaor & Draind[T993l) and use the escape probability formalism 
(Eq. As parameters at 350 and 451 days we find temper- 
atures of 700 and 600 K, and optical depths at 10 fim of 2.1 
and 3.2, respectively, to be good choices for reproducing the ob- 
served spectrum. These values are similar to the ones derived by 
K09. For the circumstellar dust component, we assume a black- 
body SED, and obtain a good fit at 350 days with T = 50 K and 
L = 3 X 10^^** ergs-i. 

We also add a component of vibrational-rotational SiO lines 
from the fundamental band at 7 - 9 fim. These are calculated 
from a NLTE solution of the SiO molecule, including the first 
100 rotational level s of the v - to v - 3 vi brati onal levels, with 
collision rates from lVarambhia et al.l (l2009h and Liu & Dalgarno' 
d 19941) . and radiative transition rates from iDrira et al.i (il997i) . 
The parameters for this component are the temperature, the elec- 
tron density, and the total luminosity. For the electron density we 



at 45 1 days reproduce the observed spectra in the 7-9 ftm range 
well. 

Figs. [To] shows the Spitzer spectrum at 350 days (red line), 
together with the 15 M© model spectrum including the ejecta 
dust, circumstellar dust, and SiO components described above 
(black line). It is clear that above 8 fim, dust emission domi- 
nates the atomic emission. At shorter wavelengths, however, the 
atomic emission is about as strong as the dust emission. It is also 
clear that CO 4.55 fim fundamental band emission is present, 
as there is not enough emission by dust and atoms in the 4-5 
fim range to account for the photometric 4.5 fim observations. 
To reach the observed flux, we need Leo = 6.0 x 10^"* erg s"' at 
350 days, and for 295 and 451 days we need Leo = 4.1 x 10^*^ 
ergs"' and Leo = 2.1x10^^ erg s"', respectively. That significant 
amounts of CO have formed is also evidenced by the presence 
of the first overtone band at 2.3 fim (Fig. |9]l. The presence of 
both CO and SiO emission supports our assumption of molecu- 
lar cooling occurring in the O/Si/S and O/C zones (see also Sec. 
|4AB. 

The atomic lines we identify using our model spectra are 
labeled in Fig.[TO| Most lines do not differ significantly between 
the models, as they mainly originate in the Fe/He or H zones. 
The only line sensitive to the core mass is [Neii] Al2.d,l fim. This 
line, arising in our models mainly from the He/C zone, but also 
from the O/Ne/Mg zone, is too strong in all models, suggesting 
that the neon content in this zone is too high, that we have mixed 
the He/C zone too strongly into the radioactive core, or that our 
ionization solution overproduces singly ionized neon. 

In our models, the [Arii] /16.985 fim line is not emitted by 
synthesized argon, but by primordial argon in the He zone. Here, 
argon is singly ionized and is responsible for about 1% of the 
cooling. The synthesized argon is located mainly in the Si/S 
zone, but is here only ionized to ~ 10^^-10"^, and provides only 
~ 10 of the cooling, which is dominated by [S i] and [Can]. 



4.4. Determination of the iron zone filling factor 

As demonstrated for SN 1987 A, the infrared iron-group 
lines are useful as diagnostics of the physical conditions in 
the ^^Ni "bubble", r esulting from radioactive ^^Ni heating 
(iMoselev et all 119891: iHaas et al.1 119901: ISpvromilio & GrahamI 
ll992l:lRoche et al.ll993HLi et al.lll993ir A complication is, how- 
ever, that there are varying contributions to these lines from pri- 
mordial iron from the H and He zones (KF98b, M12). In addi- 
tion, the studies mentioned make several assumptions about the 
ionization balance and optical depths. More reliable results re- 
quire calculation of these quantities. 

Because both the Fe/He and H/He zones are similar in the 
various models, the iron-group lines differ little between them, 
and we consider only the 15 M© model. The electron density 
in the Fe/He zone at 350 days is of order 10** cm"^ in the 
model, much higher than the critical densities for forbidden fine- 



take «e = 7 X 10 cm"^ at 350 days and «e = 2.6 x 10 cm"-* at 
451 days, as given by the models for the O/Si/S zone. Because 
of the limited S/N and resolution of the spectra, the shape of 
the P and R branches of the v = 1 to v = bands are not well 
defined. We therefore use a temperature similar to that in SN 
1987A of 1500 K (Liu & Dalgarno 19_M). We assume a velocity 
equal to that of the core velocity. For the luminosity, we find that 
Lsio = 9.9xl0"ergs-i at 350 days and Lgio = 4.3 x lO^Vg s"' = V(t)nMA- 



structure lines (~10 cm ), so all lines arising from the lower 
multiplets are in LTE. We find, for all possible filling factors of 
the iron clumps, optical depths much higher than unity for the 
[Fen] All. 94 fim and [Feii] /125.99 fim lines. These lines are 
thus both in LTE and optically thick. The general expression for 
line luminosity, in the Sobolev approximation, is 



-hv 



(9) 



12 



Anders Jerkstrand et al.: The progenitor mass of the Type IIP supernova SN 2004et from late-time spectral modeling 




10 
X [^m] 



Fig. 10. Observed MIR spectrum at 350 days (red, solid line), compared to model spectra of the 15 M© model (black, solid). To 
simulate the instrument response, the model has been convolved with a Gaussian of FWHM = A/R, where the resolving power R 
was taken to be 60. Also shown are the dust component (pink, dashed), the observed photometry at 3.6 and 4.5 j^m (red stars), and 
the corresponding model photometry for the model (black circles). The observed photometry is at 361 days, and has been rescaled 
with the ^''Co decay factor to 350 days (exp - 11 / 111 .5). 



where V{t) is the zone volume, «u(f) is the number density of the 
upper level, A is the transition probability, and Tj is the Sobolev 
optical depth. For an optically thick line in LTE, Eq. ^ becomes 



295 days 



L(t) 



V{t) 



Snhv 



t ^3(g;nvmo_ 1) 
Since hv <s: kT for the MIR lines, the line luminosity is then 



(10) 



L(t) cc A-^frT{t) , 



(11) 



where / is the filling factor Since we, at a given time t, can com- 
pute the temperature T{t) for any given filling factor /, we are 
therefore in a good position to determine the latter by comparing 
model line luminosities to the observed ones. 

Figs. fTT]fT3] show the model fluxes of the Fe/Co/Ni MIR lines 
relative to the measured fluxes, as a function of /pe/He- When we 
vary fpe/He, we add/subtract the corresponding amount to fu- We 
measure the line fluxes as described in Sect. 14.3.11 using Vcut - 
+6400kms-i(= ±3crfor7; = 60), and V^ont = l-25ycut = +8000 
kms The model spectra were convolved to the Spitzer reso- 
lution (R = 60), before this operation. For each epoch we only 
plot the lines with clear detections. 

The MIR lines behave quite consistently, many of them being 
optically thick and in LTE. The filling factor-dependency is gen- 
erally weaker than Eq. ( fTTT i suggests, as there is also an (almost 
constant) component from the He/H zones. We find that we can 
fit the Fe/Ni/Co lines only if /pe/He = 0.1 - 0.2. At larger values, 
the singly ionized lines become too strong, and in addition dou- 
bly ionized lines ([Feiii] A22.9A ^m, [Com] /ill. 89 /^m, [Com] 
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Fig. 11. Model Fe/Ni/Co line fluxes relative to observed values, 
at 295 days. The 15 Mq model has been used. 



/116.39 jim) emerge in the model. These lines are not observed, 
and we therefore rule out a filling factor larger than /pe/He ~ 0.3. 
Values much smaller than /pe/He ~ 0.1 can also be ruled out, 
as the singly ionized lines are then too weak, and unobserved 
neutral lines in the optical/NIR (Fe i] A-s{a^F)-As'\-p{z'D") tran- 
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Fig. 12. Same as Fig.ITT] for 350 days. 



an independent determination of the ^^Co mass. Also the [Con] 
/(14.74 and 15.46 /im lines have significant optical depth (t 1 
at 350 days). 

The [Ni ii] /16.634 /im line arises mainly from stable ^^'Ni in 
the Fe/He zone, with some contribution from ^^Ni. We find that 
the [Ni ii] 6.634 line is well produced by the amount of ^^Ni in 
the explosion models (M(-'''^Ni)/M(5*Ni) « 0.04). 

4.4.1. Molecules 



Figs. [T4JT5] show the estimated total luminosities of CO and SiO 
(Sect. 14.3.3b . compared to the heating of the O/C and O/Si/S 
zones in the three models. If molecules are responsible for all 
of the cooling in these zones, these values should match for the 
model with the correct O/Si/S and O/C zone masses. 

For SiO, no observations of the overtone exist, and for CO, 
we have only one observation. The CO overtone luminosity at 
the other epochs is assumed to scale to the day 307 luminosity i n 
the same manner as in SN 1987A (iBouchet & Danzigeiill993h . 
It contributes less than 10% to the total luminosity at 350 and 
450 days. 
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Fig. 13. Same as Fig.[TTl for 451 days. 



sitions between 8000 - 8500 A and [Fe i] 4s^(a^D) - 4s(a^F) 
transitions at 1.36, 1.44 and 1.53 jum), are strong. A filling factor 
of /pe/He = 0.1 - 0.2 is also supported by the [Fe 11] /I7155 line, 
whose behaviour agrees well with the MIR lines, removing any 
uncertainty regarding the MIR calibration for our conclusions. 

As is seen in the figures, [Fe 11] /124.52 /im is not reproduced 
well by any filling factor. This line is, however, weak and in a 
noisy part of the spectrum. Also the [Fe n] /125.99 fim line shows 
poorer agreement compared to the shorter wavelength lines. 

Our result here is close to the /pe/He value derive d for SN 
1987A by KF98b (/pe/He ~ 0.2). The value derived bv lLi et al.l 
( Il993h (/pe/He > 0.3) is likely too large as they ignore the pri- 
mordial Fe emission. We find that for /pe/He = 0.15 and t - 350 
days, the contribution by primordial iron-group elements in the 
H and He zones to the line fluxes is 30 - 40% for the [Fe 11] lines, 
and 20% for [Ni n] /16.636 fim. The contribution by primordial 
cobalt is low (< 1%) as the solar abundance of cobalt is very 
low. The [Co 11] /110.52 /im line thus comes from the radioactive 
^^Co, but being optically thick (up to 500 days) does not allow 
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Fig. 14. Observed flux in the CO fundamental band (hollow 
squares), first overtone (diamond), and estimated total (filled 
squares), compared to the total heating of the O/C zone in the 
various models. For the estimated total emission, correction for 
the internal dust extinction has been taken into account. All val- 
ues have been normalized to the decay of 0.062 Mq of ''^Co. 

The CO emission falls in between the 12 and 15 M© mod- 
els, behaving consistently with the optical emission lines. This 
result is also important in ruling out CO as a possibility for a 
more massive progenitor to still produce weak [O i] 6300, 6364 
A emission from the O/Ne/Mg zone. The observed CO flux in 
such a case would be much stronger 

The SiO emission is about a factor of 2-3 below the total 
heating of the O/Si/S zone in all of the models. As for CO, this 
is important for constraining the possible cooling channels of the 
O/Ne/Mg zone. Note that O/Si/S zone does not grow monotoni- 
cally with progenitor mass, being larger in the 15 M0 model than 
in the 19 Mq one (Table IS . 

It is possible that SiO only does part of the cooling of 
the O/Si/S zone. However, with the large partition function of 
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Fig. 15. Same as Fig. [14] for SiO. The overtone band is unob- 
served, and is assumed to be zero. 



molecules, one would expect it to dominate the cooling. The 
determination of the SiO fundamental band flux is complicated 
by the silicate dust feature at the same wavelength, possibly ex- 
plaining the discrepancy. It is also likely that the silicate dust, 
which presumably forms in the same zone as the SiO, produces 
stronger extinction there than for the spectrum as a whole, and/or 
does some of the cooling. 



5. Discussion 

In evaluating the robustness of linking the nebular spectra to 
a progenitor mass, there are several aspects to consider Our 
models test specific stellar evolution/explosion models (in our 
case the models by WH07), and the zone masses and com- 
positions in these may vary depending on the assumptions 



Schaller et alJl992HWooslev & Weaveill995l 


Thielemann et alj 


19961 Wooslev et alj 2002h Hirschi et al. 2004 


). It is beyond the 



scope of this paper to test all types of models with respect to 
variation in these parameters. However, an attempt to judge the 
influence of some of them is warranted. 

The influence of metallicity on the total nucleosynthesis 
yield is quite weak, although it does affect specific elements in 
particular zones (WW95). Estimates of the metallicity in the re- 
gion of SN 2004et indicate that it is likely within a factor two 
from solar (Sect. 12.21 ). and such a variation should not have any 
strong impact on the stellar evolution or SN spectral formation. 
Although it was found by Ml 2 that emission by primordial met- 
als in the hydrogen envelope can be significant contributors to 
the nebular spectrum, this emission will change little with metal- 
licity as long as the relative amounts of metals are roughly the 
same. The fraction of cooling done by each element will then 
roughly stay the same, although a slight temperature increase 
will occur for lower metallicities due to the fewer number of 
coolants available. We therefore do not think that the assump- 
tion of solar metallicity is crucial for any of our conclusions. 

Rotation can enhance the nucleosynthesis yield for a given 
stellar mass (e.g. jiirschi et al. 2004). For rapidly rotating pro- 
genitor models, our best fit would (most likely) have been a 
lower-mass star. With regard to this parameter, we can there- 



fore say that 15 is an upper limit to the progenitor mass. 
One should note, however, that the progenitor analyses (e.g. 
Crockett et al. 2011) also assume non-rotating stars, so if one 
postulates rotation, best-fitting progenitor models also go down 
in mass (rotation boosts luminosity in the late evolutionary 
stages). 

Assumptions about the symmetry, mixing, and fragmenta- 
tion of the explosion also influence the model outputs. The ob- 
served line profiles constrain the velocity range of the material, 
however, and furthermore do not indicate any extreme asymme- 
tries in the SN. As gamma-rays scatter around the whole ejecta in 
the nebular phase, neither is it obvious that even strong asymme- 
tries will significantly change the energy deposition in the vari- 
ous zones. 

Our results here are in fair agreement with the progeni- 
tor analysis of SN 2004et, which has given Mzams = H^i 
Mq (Fraser et al. 201 2, in prep.). On the other hand, hydrody- 
namical modeling by lUtrobin & Chugaj (1200^ gave an ejecta 
mass of 24.5 M©, over a factor two larger than the ejecta mass 
of our best-fitting model here (10.9 Mq), as well as any final 
mass in standard stellar evolution models (th e maximum ejecta 
mass with standard mass-loss rates is -14 Modi asen & WooslevI 
l2009l) . and will be even lower i f strong pulsation-driven mass- 
loss occurs in the RSG phase dYoon & Cantiellol l2010l) ). The 
discrepancy between the hydrodynamical modeling and the pro- 
genitor analysis/nebular-phase spectral analysis is important to 
understand. More independent analyses in all three fields are 
needed to make progress here. A possible concern regard ing the 
hydrodynamical modeling by lUtrobin & Chugail ( |2009|) is that 
non-evolutionary progenitor model s are used. 

We also mention the results in IChevalier et al.l (l2006l) . who 
estimated a progenitor mass of ~20 Mq by linking the estimated 
mass-loss rate from radio observations to standard mass-loss 
laws and stellar evolution models. This result relies on several 
assumptions (wind velocity, shell velocity, circumstellar tem- 
perature, and dominance of free-free opacity), and in addition 
the mass- l oss la ws have standard deviations of a factor 2-3. 
iRho et alj Um% derived a lower mass-loss rate of 2 x 10 ^ Mq 
yr"', which corresponds to a Mzams = 15 M© progenitor 

Comparing with other nebular-phase analysis of SN 2004et, 
we reach a similar conclusion as S06 and MIO do by simply 
comparing the [Oi] doublet to SN 1987A, that the oxygen mass 
is somewhat lower than in SN 1987 A, and the progenitor thereby 
a bit less massive (assuming similar rotation etc.). We have here, 
however, put this result on a much firmer footing by comput- 
ing detailed models and finding consistent results from several 
other lines. One should also bear in mind the contribution by 
primordial oxygen that complicates such simple comparisons 
(M12). The comparison to SN 1987A works quite well in this 
case as these two SN both have large masses of synthesized oxy- 
gen, similar velocities, and similar ^*Ni-masses. One should note 
though, that the progenitor mass of SN 1987A can not be taken 
to be known better than in the range 16-22 M© dArnett et al.l 
il989i) . Since modeling of SN 1987A would require only small 
changes to the parameters used here, we would most likely find 
a m ass on the lower side of this range. 

Dessart et al. (2010) find, using models without any macro- 
scopic mixing, a coiTelation between the width of the oxygen 
lines and the progenitor mass (for a fixed explosion energy). For 
an explosion energy of 2.3 B, they find that the progenitor mass 
of SN 2004et should be less than 18 M©, which is in agreement 
with our result here. 

For a more general comparis on of our modeling with other 
work, we consider the results in iDessart & Hilliej (1201 ih . who 
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compute model spectra for 15 and 25 Mq progenitors, also taken 
from WH07. These authors include a full NLTE and radiative 
transfer treatment, and in this sense their work is complementary 
to ours. They do, however, assume some microscopic mixing 
to occur, which may have important effects on the cooling and 
ionization balance, and therefore also on the emerging spectrum. 
At 300 days, the difference in the flux of the [O i] doublet appears 
to differ by the same amount as the ^*Ni mass in the two models 
(factor 1.9), a result possibly in disagreement with our model, 
but since we use the same ^^Ni mass for all models we cannot 
say with certainty. At 1000 days the more massive star shows 
about a factor of 3.6 stronger [Oi] AA63QQ, 6364 emission, or a 
factor 1.9 stronger per unit ^^Co mass. 

In Ml 2, it was shown that the diversity of line emission in a 
sample of Type IIP SNe is quite small. One reason for this was 
found to be that the large hydrogen envelopes in these objects 
add significant contribution to all metal emission lines, mitigat- 
ing the impact of varying core masses. Despite this damping ef- 
fect, our modeling here suggests that stars at and above 19 Mq 
should produce Nai ^^5890, 5890, [Oi] ^^6300, 6364, [Neii] 
/(12.81 fim, [Sii] /II. 60, 1.64 fim about a factor two stronger 
than lower-mass progenitors (15 Mq and below). None of the 
SNe in the Ml 2 sample, including SN 1987 A, show such strong 
oxygen emission. Future detailed modeling of individual objects 
(making sure the "missing flux" cannot be accounted for in other 
lines such as CO), and a larger sample size, will show if these 
high-mass stars are indeed missing from the Type IIP group. 

6. Conclusions 

- We have shown that nebular-phase spectral modeling can be 
used to constrain the progenitor masses of Type IIP super- 
novae. Lines showing significant dependency on the progen- 
itor mass are [Oi] AA630Q, 6364, Nai ^^5890, 5896, Mgi] 
.14571, Mgi ^1.503 fim, Ci ^1.176 fim, Ci .11.454 //m, Sii 
Al.lOfim, [Sii] .1.11.607, 1.645 //m, and [Neii] .112.81 fim. 

- The nebular-phase optical and near-infrared emission lines 
of SN 2004et are well reproduced by model spectra of a 
A^ZAMS = 15 Mq (non-rotating) progenitor star, with an 
oxygen mass of 0.8 Mq. Spectra from 12 and 19 Mq pro- 
genitors under- and overproduce most of the nebular phase 
lines, respectively, with particular discrepancy for the high- 
mass model. This result is in fair agreement with the anal- 
ysis of pre-explosion imaging (Mzams - H^jA^o, Fraser et 
al. 2012, in prep.), but in disagreement with the progenitor 
mas s inferred from hydrodynamical modeling (Mzams >25 
M^. IUtrobin&Chu gai (2009)). 

- In the 300 - 600 day range, we find that the [O i] 6300, 6364 
doublet emits 1, 2, and 4 - 6% of the ^*'Co decay energy for 
12, 15, and 19 Mq progenitor stars. This result may be used 
for a coarse diagnosis of other Type IIP SNe with similar 
5^Ni mass (0.06 Mq). 

- By modeling the mid-infrared iron-group lines observed 
with Spitzer, we determine a filling factor / = 0.15 for the 
Fe/Co/Ni clumps in the ejecta, corresponding to a density 
pit) = 7 X 10"'"^ X (f/100 d)"^ g cm"^ This is similar to the 
filling f actor deri ved for SN 19 87A by l&jzma & Fransson 
(Il998bl) . whereas in et al.l (Il993h likely overestimate the fill- 
ing factor due to neglect of emission by primordial iron in 
the hydrogen zone. 

- By our calculation of the atomic emission lines and continua, 
we can confirm the contribution by a silicate dust feature 
to the mid-infrared spectrum, first identified by iKotak etal] 
(I2009h . as well as the presence of SiO and CO fundamental 



band emission in the spectrum. The CO luminosity matches 
the total heating of the O/C zone in a 12 - 15 Mq progenitor. 
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Appendix A: Chemical composition of explosion 
models 

Tables lAniA.2l and lA.3l show the chemical composition of the 
three e xplosion models we use (taken from IWooslev & Hegej 
(120071) '). 

Appendix B: Code updates compared to J11 

B.1. Photoexcitation/deexcitation rates 

Because the computations here are for earlier epochs than in J 11 , 
photoexcitation/deexcitation rates in the NLTE solutions have to 
be included. We therefore implemented a modified Monte Carlo 
operator to allow for partial deposition of photoexcitation energy 
along the photon trajectories. A photon packet of comoving en- 
ergy E and frequency v, coming into resonance with a transition 
Z — > M of atom k, where the upper level u < u^^^ (=the number of 
levels included in the NLTE solution of atom k), now increases 
the number of photoabsorbed packages by 

AA^i,k,i,u = -;^(l-e""'"), (B.l) 
nv 

where / is the zone index and t^^ the Sobolev optical depth. The 

packet then continues (undeflected) with energy Zse^'^K For lines 
with u > we retain the treatment described in J 11 (scatter- 
ing/fluorescence is treated explicitly in the transfer). 

The values for t/j^,^^ should be chosen large enough to reason- 
ably well sample the total number of photoexcitations from each 
level, but small enough to ensure stability and keep computation 
times limited (matrix inversion times in the statistical equilib- 
rium solutions scale as mJ^.^^ cubed). The mI^,.,^ values used here 

I lid A ' IIldA 

can be found in Table IB. 21 We have chosen to set m^^^ - 300 
for all iron-group elements (except Fe i, Fe ii and Cr ii, where we 
solve for more levels), and u^^^ = 200 for the non-iron-group 
elements. 

In the Sobolev approximation, the photoexcitation rates to be 
used in the statistical equilibrium equations are given by dCastorl 
119701) 

^lu = Blu^lu-^v ' (B-2) 

where Biu is the Einstein coefficient for photoexcitation, is 
the escape probability, and is the far blue-wing mean intensity 
(before Jy is affected by the line). The deexcitation rate is related 
by detailed balance 

^ui = — ^lu , (B-3) 

where gi and are the statistical weights. 

Overlapping lines may cause to significantly vary with the 
choice of the blue frequency. We therefore do not compute the 
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Table A.l. Chemical composition (mass fractions) of the 12 Mq model used. All values are before any decay of the radioactive 
isotopes ^*Ni, ^^Co, and '*^Ti has occurred. Mass fractions below 10"*^ are set to zero. 



Element / Zone 


Fe/He 




Si/S 




O/Si/S 




O/Ne/Mg 


o/c 




He/C 




He/N 






H 








0.57 




0.058 






































^'Co 


0.034 




1.5 X 10 


_^ 


2.1 X 10" 




2.6 X 10" 


-7 


3.8 X 10" 




4.1 X 10" 




1.4 X 10" 











44 


5.9 X 10" 




1.5 X 10" 






































H 


4.0 X 10" 


— 


1.2 X 10" 


— 


2.3 X 10" 


— 


7.0 X 10" 


-I) 


2.8 X 10" 


— 


2.0 X 10" 


— 


4.7 X 


10" 


— 


0.68 






He 


0.31 




8.9 X 10" 




5.3 X 10" 




4.5 X 10" 


-6 


2.5 X 10" 




0.73 




0.99 






0.31 






C 


3.1 X 10" 




1.6 X 10" 




2.2 X 10" 




3.9 X 10" 


-3 


0.19 




0.23 




2.1 X 


10 




1.6 X 


10 




N 


3.4 X 10" 









2.9 X 10" 




3.3 X 10" 


-5 


2.1 X 10" 




3.1 X 10" 




9.1 X 


10 




2.6 X 


10" 







2.6 X 10" 




1.1 X 10 




0.80 




0.66 




0.74 




0.020 




2.2 X 


10" 




5.8 X 


10" 




Ne 


2.8 X 10" 




5.0 X 10" 




1.0 X 10" 




0.25 




0.046 




0.013 




1.1 X 


10" 




1.3 X 


10" 




Na 


9.9 X 10" 




8.2 X 10" 




1.6 X 10" 




6.2 X 10" 


-3 


1.8 X 10" 




2.3 X 10" 




1.6 X 


10 




5.6 X 


10" 




Mg 


5.4 X 10" 




1.2 X 10" 




0.044 




0.062 




0.017 




5.0 X 10" 




7.2 X 


10 




7.3 X 


10" 




Al 


6.1 X 10" 




1.9 X 10" 




3.2 X 10" 




4.0 X 10" 


-3 


1.3 X 10 




5.8 X 10" 




7.3 X 


10" 




6.6 X 


10" 




Si 


2.4 X 10" 




0.41 




0.12 




4.3 X 10" 


-3 


9.2 X 10" 




9.2 X 10" 




8.2 X 


10" 




8.2 X 


10" 




S 


2.3 X 10" 




0.39 




0.020 




3.0 X 10" 


-4 


2.2 X 10" 




3.1 X 10" 




4.2 X 


10 




4.2 X 


10" 




Ar 


2.7 X 10" 




0.054 




5.7 X 10" 




8.5 X 10" 


-5 


8.9 X 10" 




8.2 X 10" 




1.1 X 


10" 




1.1 X 


10" 




Ca 


3.1 X 10" 


-3 


0.034 




2.5 X 10" 


-5 


3.9 X 10" 


-5 


2.5 X 10" 


-5 


3.5 X 10" 


-5 


7.4 X 


10" 


-5 


7.4 X 


10" 


-5 


Sc 


2.1 X 10" 


-7 


2.3 X 10" 


-7 


6.7 X 10" 


-8 


8.8 X 10" 


-7 


9.1 X 10" 


-7 


1.0 X 10" 


-6 


5.6 X 


10 


-8 


4.5 X 


10" 


-8 


Ti(stable) 


1.9 X 10" 


-3 


5.3 X 10" 


-4 


7.5 X 10" 


-6 


5.3 X 10" 


-6 


8.0 X 10" 


-6 


3.2 X 10" 


-6 


3.4 X 


10 


-6 


3.4 X 


10" 


-6 


V 


6.1 X 10" 


-5 


2.2 X 10" 


-4 


2.4 X 10" 


-6 


3.7 X 10" 


-7 


3.7 X 10" 


-7 


6.0 X 10" 


-7 


4.5 X 


10 


-7 


4.3 X 


10" 


-7 


Cr 


2.9 X 10" 


-3 


6.9 X 10" 


-3 


9.1 X 10" 


-6 


1.6 X 10" 


-5 


1.2 X 10" 


-5 


1.5 X 10" 


-5 


2.0 X 


10 


-5 


2.0 X 


10" 


-5 


Mn 


1.9 X 10" 


-6 


2.9 X 10" 


-4 


1.2 X 10" 


-6 


7.6 X 10" 


-6 


1.9 X 10" 


-6 


3.9 X 10" 


-6 


1.6 X 


10 


-5 


1.5 X 


10" 


-5 


Fe (stable) 


7.0 X 10" 


-4 


0.045 




3.9 X 10" 


-4 


9.2 X 10" 


-4 


5.1 X 10" 


-4 


1.1 X 10" 


-3 


1.4 X 


10" 


-3 


1.4 X 


10" 


-3 


Co (stable) 


2.2 X 10" 


-8 


1.9 X 10" 


-9 


1.2 X 10" 


-4 


1.7 X 10" 


-4 


1.8 X 10" 


-4 


2.5 X 10" 


-4 


4.3 X 


10" 


-6 


4.0 X 


10" 


-6 


Ni(stable) 


0.029 




2.5 X 10 


-3 


9.5 X 10 


-4 


3.9 X 10" 


-4 


7.2 X 10" 


-4 


1.4 X 10" 


-4 


8.2 X 


10" 


-5 


8.2 X 


10" 


-5 



Table A.2. Same as Table [Al] for the 15 Mq model. 



Zone 


Fe/He 




Si/S 




O/Si/S 




O/Ne/Mg 


O/C 




He/C 




He/N 






H 






^"Ni + ^"Co 


0.66 




0.057 




1.9 X 10 


-7 


3.5 X 10" 


-7 


1.1 X 10" 


-8 


7.5 X 10" 


-8 


2.0 X 10" 


-8 









5^Co 


0.033 




1.4 X 10" 


-3 


3.1 X 10 


-6 


2.6 X 10" 


-7 


1.3 X 10" 


-8 


6.6 X 10" 


-9 
















44 


5.1 X 10" 


-4 


1.5 X 10" 


-5 


1.1 X 10 


-9 































H 


2.5 X 10" 


-6 


4.0 X 10" 


-7 


7.4 X 10" 


-9 


1.6 X 10 


-9 












4.5 X 


10 


-8 


0.65 






He 


0.23 




9.5 X 10" 


-6 


4.0 X 10" 


-6 


3.7 X 10" 


-6 


1.6 X 10" 


-4 


0.94 




0.99 






0.33 






C 


1.4 X 10" 


-6 


1.3 X 10" 


-6 


3.9 X 10" 


-4 


8.5 X 10" 


-3 


0.20 




0.039 




2.4 X 


10 


-4 


1.5 X 


10" 


-3 


N 


2.3 X 10" 


-6 







2.9 X 10" 


-5 


6.9 X 10 


-5 


1.3 X 10" 


-5 


2.7 X 10" 


-3 


9.1 X 


10 


-3 


3.0 X 


10- 


-3 





1.7 X 10" 


-5 


1.0 X 10 


-5 


0.81 




0.68 




0.73 




5.6 X 10 


-3 


1.8 X 


10" 


-4 


5.4 X 


10- 


-3 


Ne 


1.9 X 10 


-5 


7.3 X 10" 


-6 


1.4 X 10" 


-4 


0.23 




0.050 




6.8 X 10" 


-3 


1.1 X 


10" 


-3 


1.2 X 


10- 


-3 


Na 


8.6 X 10" 


-7 


7.9 X 10" 


-7 


2.1 X 10 


-5 


5.2 X 10" 


-3 


1.9 X 10 


-4 


1.8 X 10" 


-4 


1.8 X 


10 


-4 


6.1 X 


10" 


-5 


Mg 


4.0 X 10" 


-5 


1.4 X 10" 


-4 


0.045 




0.062 




0.016 




7.3 X 10" 


-4 


7.0 X 


10 


-4 


7.3 X 


10" 


-4 


Al 


1.0 X 10 


-5 


2.0 X 10" 


-4 


4.4 X 10" 


-3 


4.0 X 10" 


-3 


1.2 X 10" 


-4 


7.2 X 10" 


-5 


9.0 X 


10" 


-5 


6.6 X 


10" 


-5 


Si 


2.4 X 10" 


-4 


0.41 




0.12 




4.8 X 10" 


-3 


9.4 X 10" 


-4 


8.2 X 10" 


-4 


8.2 X 


10" 


-4 


8.2 X 


10" 


-4 


S 


2.2 X 10" 


-4 


0.39 




0.019 




2.9 X 10" 


-4 


2.2 X 10" 


-4 


4.2 X 10" 


-4 


4.2 X 


10" 


-4 


4.2 X 


10" 


-4 


Ar 


2.3 X 10" 


-4 


0.055 




5.9 X 10" 


-4 


8.4 X 10" 


-5 


8.6 X 10 


-5 


1.1 X 10 


-4 


1.1 X 


10 


-4 


1.1 X 


10- 


-4 


Ca 


2.7 X 10" 


-3 


0.035 




2.8 X 10" 


-5 


3.6 X 10" 


-5 


2.6 X 10" 


-5 


7.3 X 10" 


-5 


7.4 X 


10" 


-5 


7.4 X 


10- 


-5 


Sc 


2.2 X 10" 


-7 


2.2 X 10" 


-7 


1.2 X 10" 


-7 


1.2 X 10" 


-6 


1.6 X 10" 


-6 


7.3 X 10" 


-8 


4.5 X 


10" 


-8 


4.5 X 


10- 


-8 


Ti (stable) 


1.6 X 10 


-3 


5.6 X 10" 


-4 


8.1 X 10 


-6 


5.6 X 10" 


-6 


7.0 X 10" 


-6 


3.4 X 10" 


-6 


3.4 X 


10 


-6 


3.4 X 


10" 


-6 


V 


3.3 X 10" 


-5 


1.5 X 10" 


-4 


3.0 X 10" 


-6 


4.1 X 10" 


-7 


4.5 X 10" 


-7 


4.7 X 10" 


-7 


4.3 X 


10 


-7 


4.3 X 


10" 


-7 


Cr 


2.3 X 10" 


-3 


6.6 X 10" 


-3 


8.2 X 10" 


-6 


1.5 X 10 


-5 


1.2 X 10" 


-5 


2.0 X 10" 


-5 


2.0 X 


10" 


-5 


2.0 X 


10" 


-5 


Mn 


1.7 X 10" 


-6 


2.8 X 10" 


-4 


1.7 X 10 


-6 


6.7 X 10" 


-6 


2.2 X 10" 


-6 


1.7 X 10" 


-5 


1.5 X 


10" 


-5 


1.5 X 


10" 


-5 


Fe (stable) 


7.6 X 10" 


-4 


0.042 




3.4 X 10" 


-4 


8.5 X 10" 


-4 


5.5 X 10" 


-4 


1.4 X 10" 


-3 


1.4 X 


10 


-3 


1.4 X 


10- 


-3 


Co (stable) 


2.3 X 10" 


-8 


2.0 X 10" 


-9 


1.9 X 10 


-4 


1.8 X 10" 


-4 


2.0 X 10" 


-4 


4.8 X 10" 


-6 


4.0 X 


10 


-6 


4.0 X 


10" 


-6 


Ni (stable) 


0.028 




2.3 X 10" 


-3 


8.6 X 10 


-4 


4.4 X 10" 


-4 


6.7 X 10" 


-4 


8.2 X 10" 


-5 


8.2 X 


10 


-5 


8.2 X 


10" 


-5 



photoexcitation rates from this equation, but by adding the con- 
tributions from Eq. (IB. lb . Since stimulated emission is treated as 
negative absorption in the Sobolev optical depth, Eq. ( IB. II ) gives 
only the net number of photons absorbed in the line; absorptions 
minus stimulated emissions. The actual number of absorptions 
occurring is therefore higher by afactor(l -gi«u/§uni) where 



«i and «u are the populations of levels / and u, so the estimator 
for the photoexcitation rate is 

, (B.4) 

ni,k,i \ ^k,uni,k,i/ 

where A^i,k,i,u is the total number of packets absorbed in the tran- 
sition in the Monte Carlo simulation. 
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Table A.3. Same as TablelAl] for the 19 Mq model. 



Zone 


Fe/He 








U/ai/J) 




U/iNe/Mj 




u/c 




He/C 




He/N 






TT 
H 






JNi+ Co 


0.76 




0.029 




5.8 X 10" 


-7 


1.3 X 10" 


-7 


9.4 X 10" 


-ij 


5.0 X 10" 


-8 


1.6 X 10 


-8 









57 

Co 


0.033 




1.0 X 10" 


-3 


1.1 X 10 


-5 


9.8 X 10 


-8 


1.9 X 10" 


-8 


3.7 X 10" 


-9 


















2.5 X 10" 


-4 


1.1 X 10" 


-5 


4.8 X 10" 


-6 































H 


2.5 X 10" 





4.1 X 10" 


— 


8.6 X 10" 


— 

















3.7 X 


10" 





0.62 






He 


0.13 




8.5 X 10" 


-6 


4.3 X 10" 


-6 


2.0 X 10" 




0.078 




0.96 




0.99 






0.37 






C 


3.2 X 10" 


-7 


1.3 X 10" 


-6 


8.5 X 10" 


-5 


0.031 




0.28 




0.023 




2.6 X 


10" 


-4 


1.4 X 


10" 


_3 


N 


1.4 X 10" 


-f, 







1.8 X 10" 


-5 


3.2 X 10" 


-5 


4.8 X 10" 


-f, 


6.4 X 10" 


-4 


9.0 X 


10" 


-3 


3.7 X 


10 


-3 


O 


9.2 X 10 


-6 


1.3 X 10" 


-5 


0.45 




0.70 




0.62 




3.3 X 10" 


-3 


1.7 X 


10" 


-4 


4.9 X 


10" 


.3 


Ne 


9.6 X 10 


-6 


1.6 X 10" 


-6 


1.4 X 10" 


-4 


0.22 




0.017 




0.011 




1.1 X 


10" 


-3 


1.2 X 


10" 


-3 


Na 


6.2 X 10" 


-7 


1.2 X 10" 


-6 


1.2 X 10" 


-6 


7.0 X 10" 


-3 


1.9 X 10" 


-4 


1.8 X 10" 


-4 


1.8 X 


10" 


-4 


7.2 X 


10 


-5 


Mg 


8.9 X 10" 


-6 


1.7 X 10" 


-4 


7.2 X 10" 


-4 


0.042 




4.9 X 10" 


-3 


7.2 X 10" 


-4 


7.0 X 


10" 


-4 


7.3 X 


10 


-4 


Al 


2.1 X 10 


-5 


2.6 X 10" 


-4 


1.9 X 10" 


-4 


4.0 X 10" 


-3 


6.7 X 10" 


-5 


7.5 X 10" 


,5 


9.8 X 


10" 


-5 


6.7 X 


10" 


.5 


Si 


1.5 X 10 


-4 


0.45 




0.27 




3.1 X 10" 


-3 


8.9 X 10" 


-4 


8.2 X 10" 


-4 


8.2 X 


10" 


-4 


8.2 X 


10 


-4 


S 


1.1 X 10 


-4 


0.39 




0.21 




2.6 X 10" 


-4 


3.2 X 10" 


-4 


4.2 X 10" 


-4 


4.2 X 


10" 


-4 


4.2 X 


10 


-4 


Ar 


1.0 X 10" 


-4 


0.049 




0.043 




8.4 X 10" 


-5 


8.9 X 10" 


-5 


1.1 X 10 


-4 


1.1 X 


10" 


-4 


1.1 X 


10" 


-4 


Ca 


1.3 X 10" 


-3 


0.029 




0.013 




3.1 X 10" 


-5 


5.0 X 10" 


-5 


7.3 X 10" 


-5 


7.4 X 


10" 


-5 


7.4 X 


10 


.-5 


oc 


Z.D X iU 


-7 


Z. 1 X iU 


-7 


1 .Z X iU 


-6 


1 A V 1 
i .0 X iU 


-6 


J.J X iU 


-7 


/.OX iU 


-8 


4.5 X 


10" 


-8 


4. J X 


1 r\' 

iU 


-8 


Ti (stable) 


8.1 X 10" 


-4 


3.7 X 10" 


-4 


9.3 X 10" 


-5 


6.2 X 10" 


-6 


4.4 X 10" 


-6 


3.4 X 10" 


-6 


3.4 X 


10" 


-6 


3.4 X 


10" 


-6 


V 


2.1 X 10 


-5 


1.5 X 10" 


-4 


3.7 X 10" 


-6 


4.1 X 10" 


-7 


3.4 X 10" 


-7 


4.8 X 10" 


-7 


4.3 X 


10" 


-7 


4.3 X 


10 


-7 


Cr 


1.3 X 10" 


-3 


4.5 X 10" 


-3 


4.5 X 10" 


-5 


1.3 X 10" 


-5 


1.7 X 10" 


-5 


2.0 X 10" 


-5 


2.0 X 


10" 


-5 


2.0 X 


10 


-5 


Mn 


2.7 X 10" 


-6 


2.7 X 10" 


-4 


2.2 X 10" 


-6 


3.8 X 10" 


-6 


9.4 X 10" 


-6 


1.7 X 10" 


-5 


1.5 X 


10" 


-5 


1.5 X 


10 


-5 


Fe (stable) 


1.1 X 10" 


-3 


0.050 




4.5 X 10" 


-4 


6.9 X 10" 


-4 


1.1 X 10" 


-3 


1.4 X 10 


-3 


1.4 X 


10" 


-3 


1.4 X 


10" 


-3 


Co (stable) 


3.6 X 10" 


-8 


4.4 X 10" 


-9 


1.4 X 10" 


-6 


1.9 X 10" 


-4 


1.3 X 10" 


-4 


4.4 X 10" 


-6 


4.0 X 


10" 


-6 


4.0 X 


10" 


-6 


Ni (stable) 


0.038 




2.4 X 10" 


-3 


1.1 X 10" 


-3 


5.6 X 10- 


-4 


2.3 X 10" 


-4 


8.2 X 10" 


-5 


8.2 X 


10" 


-5 


8.2 X 


10 


-5 



For lines with t^^^ < 0, we allow lasering to occur in the trans- 
fer by using the same formalism as for normal positive optical 
depths (Eq. dB.lb ). In these cases we compute the deexcitation 
rate as 



|A^i.k.i,ul 



i,k,u,l 



^k,u«i.k.l 



ni.k.u \ gk,ini,k,u I 

and then use Eq. ( IB.3b to get the excitation rate. 



B.2. Charge transfer 

We now generalize the treatment of the O n -H Ci — »0 i(2/?('D)) 
-H Cn reaction in Jll, by assuming that all charge transfer re- 
actions occur to the state (either in the recombining or in the 
ionizing element) that gives the smallest energy defect. 

B.3. Photoionization rates and photoelectric heating rates 

We have also modified the estimators for the photoionization 
rates and the photoelectric heating rates. Whereas the previous 
treatment did not allow splitting of packets upon continuum in- 
teraction, a photon packet passing a region of optical depth r, of 
which photoionization by atom k, level /, contributes r^:', now 
adds a contribution 



m.u. 



i,k,l 



^k,l 
T . p 

T llV 



(B.6) 



to the number of photons Nix,\ absorbed in the corresponding 
continuum. The heating rate is increased by 



k,l 



(B.7) 



where Vi is the volume of the zone, and xy,\ is the ionization 
potential. 



B.4. Electron scattering 

We implement electron scattering in the following way. The to- 
tal continuum optical depth to the next point of interest (a zone 
boundary, a line, or a photoionization threshold) is 



(B.5) ATcont = Arpi + Ard + Att , 



(B.8) 



where Arpi is the photoionization optical depth, Ard is the dust 
optical depth, and Att is the Thomson optical depth. The frac- 
tion of photons absorbed in photoionization continua is 



^ * cont 

the fraction absorbed by dust is 

ATd 



/c 



At,, 



the fraction that Thomson scatter is 

ATcont ^ ' 

and the fraction that do not interact is 

/pass = e"*'"™"' = 1 - /pi - /d - /j . 



(B.9) 



(B.IO) 



(B.ll) 



(B.12) 



We let the fractions /pi and /d be absorbed in photoionization 
continua and by dust, respectively. We then let the remaining 
fraction (/t + /pass) Thomson scatter with probability 



Pi = 



/t 



/t + /pass 



(B.13) 



and continue through with probability 1 - pjAf Thomson scat- 
tering occurs, we draw a position for the scattering from 



At = -ln[l +z(e"^^™°' - 1)] , 



(B.14) 
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and a new direction cosine /i from 



z - 



where z is a uniform random number between and 1, and p(iu') 
is the Rayleigh phase function 



pO"') = t(i+/^") 



(B.16) 



We treat the scattering as coherent in the fluid frame, i.e. 
we ignore the small energy shifts occurring due to the thermal 
motions of the electrons, as well as the Compton shift. 



« = 3 and O i 3d^D is assumed, the branching probability for O i 
fluorescence ifl 



<B'5) A„ =0.043^ 



(B.19) 



where n(0 1) and n(H i) are the number densities of O i and H i. 
At solar metallicity n(0)/n(H) = 4.9 x 10"^^ ^l^^ = 2.1 x 10"^ 
which is the value we use. 

B.6. Photoionization cross sections 

For excited states with unknown photoioni zation cross sections, 
we u se the hydrogenic approximation dRvbicki & LightmanI 
119791) 



B.5. Line overlap 

When the density is high, the line escape probabilities may be 
enhanced from their Sobolev values by overlapping continua 
and lines, so that 



/r« +if 



(B.17) 



We take the ff^ term into account as described in Jll. In this 
paper, we also take the jS^ term into account for hya and LyyS, 
which are important for the spectral formation in Type II SNe. 

Lye reaches optical depths as high as = 10'" - 10" in 
the early nebular phase, and may branch into several iron-group 
lines around 1215.67 A before escaping the resonance layer. 
Preliminary simulations (which we hope to expand on in a future 
publication) show that at typical conditions in the early nebular 
phase, the branching probability is 



(B.18) 



Here we use 0^^^ = 5 x 10 ^. We assume that the absorbing 
lines are the ones which, at 300 days in the core H region, are 
optically thick and within 100 Doppler widths (on the red side) 
of the Lya rest wavelength, with allocation to these transitions 
in equal proportions. The lines, which fall in the 1215.67 - 1220 
A range, are listed in Table IB. ll 



Table B.l. Transitions into which (part of the) Lyce flux is allo- 
cated. 



Fell transitions: 
1217.85 As{a*D-ji2) - (5£))5p9/2 
1217.15 As{a^Dii2) - Api^PI/^) 
1215.98 As{a^Diii) - 4pC*5°/2) 
1215.85 45(a'*£)5/2) - 5p(''£>°/2) 
1216.24 As(a*Di,2) - 5p{*Pl^,) 
1216.27 As{a*Dii2) - Api^P^i^) 
1216.52 4i(a*Z)i/2) - 5p(*D1i^) 
1216.85 AsQj^F^ii) - ^P°Ch1^^^) 
1217.62 As(b*FTi2) - ^P-'C^H"^!^) 
1218.09 4s{a^D3,2} - Spi^'F^,^) 
1218.19 4i-(fo''F5/2)- 



SI2' 

'P°{'H-i2) 



1218.70 As{ci*Dii2) - SpCP^^,) 



CTki =7.9x 10" 



1219.41 'idila^G^a) - 4p('//°„) 
1219.80 As(a*Hni2) 



gk,i«k,iZi^ cm- , 



(B.20) 



where v°j is the ionization threshold, is the Gaunt factor 
(which we set to unity for all levels), is the effective nuclear 
charge (number of protons minus number of screening elec- 
trons), and «k,i is the effective principal quantum number 



«k,l = I 1 



-1/2 



(B.21) 



where /Skj is the level energy and Xiifi is the ground state ioniza- 
tion potential. 

We use a threshold for how many levels to compute pho- 
toionization cross sections for. For this paper this limit was set 
to 50 levels. Tests showed that contributions by levels above this 
were small for all epochs modeled. 

6. 7. Free-bound emissivities 

For hydrogen, we use detailed emissivity functions obtained by 
applying the Milne relation s to the p hotoionization cross sec- 
tions, which we obtain f rom lVerner et al. (1996i) for n - I, and 
iBethe & Salpeteil ( ll957h for n - 2. For n - 3 and higher we use 
the hydrogenic approximation. 

For all other atoms k we use a simplified prescription, where 
recombination to level I gives an emissivity of 



4nj[f = «en+ak,i(r)Kk,i(r)-^r + 



(B.22) 



9/: 

13/2;- ^lip) 
Cm transitions: 
1216.31 As{a^Dii2) - 5p(^Pi,2) 
1217.31 Asia^Dyi) - spix^DjpJ 
1218.63 4s(a^Dy2) - Sp^PipJ 
1217.14 4sia^Dg/2) - spix^Dgp) 
1217.85 4j(a«£)3/2) - Spi'-PspJ 
1218.91 45(fl«£»,/2)-5p(«/'3/2) 
1219.56 41VD5/2) - SpCPsp) 
1219.96 4iVZ)3/2) - 5p(^P3pJ 



The Oi /11025.76 2//(V) - Sd^D") transition lies within 
the Doppler core of Ly/3. If the Oi /11.129 pm 3pCP) - 3dCD°) 
transition is optically thin, and complete /-mixing in both Hi 



IS emitted over a box-shaped profile between 0.5 - 1 .5 kT. Here 
(^k.\{T) is the specific recombination coefficient, and ^k,i(r) is 
taken as 0.55 for all atoms, levels , and tempera tures, the value 
for hydrogen at a few thousand K (ISpitzeijll948h . 

6.8. Coliision strengtlis 

For transitions with u nknown collision stren gths we use van 
Regemorter's formula jvan Regemorteij[l962l) for allowed lines 
(taken as /osc > lO"-'), and 

Oiu = 0.004gig, , (B.23) 

for forbidden lines (lAxelrodlll980l) . 

' 5/9 of Oi can absorb the line, there is a ~30% chance that fluores- 
cence occurs in O I, the Ag^ / gi values are 1.3 x 10* for the O i line and 
5.0 X 10' forLyjS^ /> = 5/9 -0.3 • 1.3 x 10^/5.0 x lO** = 0.043. 
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B.9. Model atoms 

Table IBT2] lists the model atoms used. Compared to Jll, we 
now include more levels in the NLTE solutions, and have also 
adde d Fe iii and Co nfl We hav e added collision strengths for 
Nin dCassidv et al.l 120101 1201 ll we use the values at 5000 K). 
The sources for the rest of the atomic data are given in Jll. 

Table B.2. Model atoms used. The Levels column shows the 
number of levels (including spin-orbit splitting) included in the 
NLTE solutions (Mmax), followed by the number included in the 
radiative transfer. The column to the right identifies the corre- 
sponding levels. See Jll for references for the atomic data. 



Element 


Levels 


Highest term 


Hi 


210/465 


n=20 / n=30 


He I 


29/29 




Ci 


90/90 


SdCD") 


Cn 


77/77 




Ni 


68 / 68 


5s(-P) 


Nil 


148 / 148 


SpiD") 


Oi 


135 / 135 




On 


16/16 


2s.2p*CS) 


Nei 
Nen 


3/3 
3/3 


ip^Cs) 


Nai 


66/66 




Mgi 


200 / 371 


37dCD)/S0dCD) 


Mg n 


35/35 


IpCP") 


All 


86/86 


2\s{^S) 


Alii 


154 / 154 


1 w.,'! J/S T\0\ 

ipid{ D ) 


Si I 


200 / 494 


\Qs{^ P)l2\s{^ P) 


Sin 


77 / 77 


9p{-P°) 


Si 
Sn 
An 


125 / 125 

5/5 
3/3 


3p^(^P°) 
4.s(-[3/2]") 


Am 


2/2 


Q '1/2 no\ 

ip ( P ) 


Cai 


198/198 


\6dCD) 


Can 


69/69 


Idd^D) 


Sci 


254 / 254 


lp(-P") 


Sen 


165 / 165 




Til 
Tin 


300 / 394 
212/212 


5pCF)/4;r(e'P) 
5di*H) 


Vi 
Vn 
Cri 


300 / 502 
300 / 323 
300 / 392 


4p(t'*G")/4p\r-H'') 

Fii24fein/2r)/Fgi24fe[n/2]") 

Ap{t^P°)l2lp(e^F) 


Crn 


364"/ 500 


4s4p(x''D°)/4fCH°) 


Mm 


293 / 293 


4p(z^n 


Mnn 
Fei 


300 / 509 
496 / 496 


4s4p(w'D")/9hCH'') 
4s4p(pH'') 


Fen 


533*/ 578 


4s4p(*H'')/4d{*D) 


Fein 


300 / 600 


*P0F'')/5pCF") 


Co I 


300/317 


4d(e^H)l5s{'iT) 


Con 


108/108 


Id^H) 


Com 


300 / 306 


4d(e*F)l5s(fD) 


Nil 


136/136 


4rf(-[l/2]) 


Nin 


300 / 500 


5p(,*S'')lbgem) 



Notes. Lya pumping occurs up to level 364. Lyo- pumping occurs 
up to level 533. 



Data from http : //kurucz .harvard, edu/atoms .html 



